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ABSTRACT 


A  reverberation  model  based  on  the  parabolic  approximation  is  developed  that 
includes  sediment  interface  and  volume  perturbations.  A  multiple  forward/single 
backscatter  approximation  is  made,  and  the  structure  of  the  solution  is  found  to  depend 
on  the  two-way  propagation  with  a  scattering  strength  scaling  dependent  on  the  local 
properties  of  the  perturbation.  The  model  is  implemented  for  continuous  wave  (CW) 
signals  to  predict  mean  reverberation  pressure  levels  and  for  broadband  pulse  signals  to 
generate  complex  reverberation  structures  in  the  time-domain.  The  spatial  correlation  and 
statistical  properties  of  these  predicted  signals  are  then  analyzed  in  an  attempt  to  extract 
information  on  the  underlying  characteristics  of  the  perturbation.  Preliminary  analysis 
suggests  that  reverberation  due  to  the  volume  perturbations  decorrelates  more  rapidly 
over  depth  than  the  reverberation  due  to  interface  fluctuations,  although  the  differences 
appear  small.  Additionally,  the  statistical  character  of  the  reverberation  structure  due  to 
the  interface  appears  as  a  relatively  flat  spectrum,  while  the  spectrum  of  the  volume 
reverberation  tends  to  appear  colored.  Attempts  to  correlate  these  characteristics  with  the 
structure  of  the  perturbations  is  ongoing. 
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I.  INTRODUCTION 


In  1993,  the  Office  of  Naval  Research  (ONR)  sponsored  a  reverberation  program, 
the  Acoustic  Reverberation  Special  Research  Program  (ARSRP),  to  study  the  primary 
causes  and  nature  of  acoustic  reverberation  in  the  deep  ocean.  The  region  studied  was 
near  the  Mid- Atlantic  Ridge,  an  area  where  bottom  characteristics  are  dominated  by 
highly  variable  topography  and  hard,  basaltic  rock  structures  (Tucholke,  et  al,  1993a  and 
1993b).  The  dominant  feature  contributing  to  long-range  reverberation  was  then  the 
water/bottom  interface.  Furthermore,  the  general  structure  of  the  reverberation  signal  was 
found  to  correlate  quite  well  with  the  forward  propagation  (Smith,  et  al.,  1996,  and 
Makris,  et  al.,  1994),  i.e.  high  reverberation  returns  were  found  to  occur  where  energy 
interacted  strongly  with  the  bottom. 

Previous  analysis  with  a  broadband  parabolic  equation  (PE)  model  has  shown 
good  agreement  with  predicted  arrival  time  structures  due  to  deep-ocean  bottom  interface 
reverberation  (Smith  and  Cushman,  1997).  Application  of  this  approach  for  both  single 
frequency  (CW)  and  broadband  pulse  propagation  calculations  showed  good  agreement 
with  measured  data.  In  addition,  the  broadband  reverberation  predictions  indicated  that 
some  features  of  the  reverberation  signal  required  the  discrimination  of  multipath 
propagation  in  order  to  fully  resolve  the  structure,  particularly  in  CW  shadow  zones. 

A  reverberation  experiment  is  now  being  organized  by  ONR  in  the  littoral  regions 
of  the  East  and  South  China  Seas  called  ASIAEX.  In  contrast  to  the  previous  deep  ocean 
regions  studied,  this  area  has  a  much  smoother  and  softer  water/bottom  interface.  It  is 
unclear  at  this  point  whether  the  dominant  mechanism  for  bottom  reverberation  will  be 
interface  roughness  or  volume  inhomogeneities  in  the  sediment.  One  may  expect  that  the 
measured  reverberation  will  contain  statistical  information  related  to  the  character  of  the 
scattering  mechanism.  However,  the  reverberation  signal  from  the  sediment  volume  may 
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be  affected  by  forward  scatter  from  the  interface  and  vice-versa,  thereby  complicating  the 
distinction  of  the  dominant  mechanism. 

The  purpose  of  the  present  work  is  to  establish  a  model  for  predicting  the 
influence  of  propagation  on  both  interface  and  volume  reverberation  in  a  shallow  water 
environment.  This  model  is  similar  to  the  previous  PEREV  model  (Tappert  and  Ryan, 
1989),  although  it  formally  treats  the  scattering  influence  in  the  context  of  the  PE 
approximation.  It  is  further  expanded  to  treat  volume  scattering. 

Starting  with  based  on  the  Monterey-Miami  Parabolic  Equation  (MMPE) 
propagation  model  developed  by  Smith  and  Tappert  (1994)  and  Smith  (2000),  bottom 
interface  and  sediment  sound  speed  perturbations  were  incorporated  and  broadband  and 
CW  analysis  conducted.  The  theoretical  basis  of  reverberation  due  to  the  bottom  interface 
and  bottom  volume  is  developed  in  Chapter  II.  Specifically,  we  were  able  to  use  a  CW 
one-way  propagation  model  to  show  that  the  reverberant  field  is  the  coupling  of  two,  one¬ 
way  propagations.  The  theoretical  basis  for  reverberation  was  expanded  to  time-domain 
analysis  for  both  the  interface  and  volume  in  order  to  provide  simulation  for  pulsed 
propagation.  In  this  method,  it  was  shown  that  the  two,  one-way  propagations  from  the 
source  to  the  scattering  patch  and  from  the  receiver  to  the  scattering  patch  would  give  the 
reverberation  field  by  convolution  of  both  fields  in  time.  Chapter  II  also  details  the 
models  used  for  generating  realizations  based  on  statistically  controlled  random 
fluctuations  for  the  interface  and  volume,  and  offer  a  practical  implementation  that  was 
incorporated  into  the  MMPE  model.  This  new  model  is  known  as  MMPEREV  (see 
Chapter  HI). 

Secondly,  after  incorporating  random  fluctuations  of  pertinent  parameters  into  the 
MMPE  model,  post-processing  analysis  of  these  results  in  MATLAB  was  performed. 
Chapter  IV  describes  the  post-processing  and  signal  processing  implemented  in 
MATLAB  and  provides  a  visual  presentation  of  the  effects  of  reverberation.  A  single  set 
of  environmental  parameters  was  used  in  all  the  simulations  to  afford  a  controlled 
computation  of  the  effects  of  reverberation,  providing  a  qualitative  description  of  the 
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analysis.  Comparisons  between  perturbed  and  unperturbed  data  were  conducted  and 
detailed  in  Chapter  IV.  Post-processing  and  signal  analysis  included: 

•  time-domain  analysis  and  comparison  with  unperturbed  results; 

•  CW  analysis  and  comparison  with  unperturbed  results; 

•  temporal  and  range  peak  auto-  and  cross-correlations  in  depth; 

•  spectral  analyses  of  the  reverberation  and  reverberation  loss. 

This  model  will  serve  as  a  tool  for  understanding  propagation  influences  in  future 
reverberation  experiments. 
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II.  THE  REVERBERATION  PROBLEM  -  SHALLOW  WATER 
REVERBERATION  AND  SCATTERING  THEORY 


In  this  chapter,  we  present  a  theoretical  treatment  of  the  effects  of  bottom 
reverberation  and  scattering.  We  acknowledge  that  in  shallow  water,  the  effects  of 
reverberation  are  more  involved  due  to  the  thin  layer  of  water  and  the  somewhat  invariant 
sound  speed  profile  in  the  water  column.  However,  perturbations  to  the  interface  and 
sediment  volume  (hereafter  known  as  just  the  ‘volume’)  exist  and  may  be  modeled 
statistically.  These  will  be  treated  in  detail.  One  element  lacking  at  this  point  is  the  study 
into  density  variabilities  in  the  volume.  Implementation  for  density  variability  will  form 
part  of  the  scope  of  follow-on  work,  which  is  discussed  in  Chapter  V. 

The  interface  and  volume  perturbations  may  be  treated  separately  because  of  the 
nature  of  their  perturbations  -  the  interface  primarily  being  perturbed  spatially  while  the 
volume  variabilities  involve  both  the  sound  speed  and  density.  In  order  to  effectively 
model  their  total  effects,  all  interface  and  volume  perturbations  must  be  jointly 
incorporated  into  the  model  at  the  same  time.  This  has  been  integrated  into  the  MMPE 
Model  developed  by  Smith,  2000  (originally  known  as  the  UMPE  Model  developed  by 
Smith  and  Tappert,  1994). 

A.  BOTTOM  INTERFACE  SCATTERING 

We  begin  with  the  standard  wave  equation,  also  known  as  the  Helmholtz  wave 
equation.  It  is  written  for  acoustic  pressure,  p,  as 

v2p+*0Vp=o  (2>1) 

with  the  following  definitions  for  wavenumber  and  index  of  refraction,  including  the 
reference  sound  speed, 

k0=O)lc0  and  n(x)  =  c0  /  c(x) .  (2-2) 
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Now,  consider  a  rough  bottom  where  the  spatial  fluctuations  may  be  defined 
statistically.  The  boundary  between  the  water  column  and  the  bottom  is  then  defined  by, 

n2  =n2w(l-H(z-zb))  +  n2bH(z-zb),  (2.3) 

where  nw  and  nb  represent  the  spatially  varying  indices  of  refraction  in  the  water  column 

and  bottom  (volume),  respectively,  and  Zb  is  the  bathymetric  depth  at  the  interface.  The 
definition  of  the  Heaviside  step  function  used  at  the  water/bottom  interface  is  given  by 


H(z)  = 


'2’ 


0, 


z  >0 

z=0  ,  and  ^-H(z)  =  S(z). 


(2.4) 


We  treat  the  influence  of  the  roughness  of  the  bottom  interface  by  introducing 
perturbations  to  the  bottom  depth  of  varying  length  scales.  This  is  akin  to  the  method  of 
small  perturbations  (MSP)  found  in  Brekhovskikh  and  Lysanov  (1990).  More 
specifically,  we  define 

Zb(r)  =  zb(r)  +  tjl(r)  +  T]s(r),  (2.5) 

where  zb  (?)  is  the  deterministic  portion  of  the  bathymetry  as  a  function  of  horizontal 
range,  F .  The  deterministic  part  of  the  bathymetry  may  also  be  seen  as  the  mean  depth  on 
a  relatively  flat  portion  .of  the  seabed.  Correspondingly,  tjt  (?)  is  the  long-wavelength 
portion  of  the  randomly  generated  roughness  comprising  a  spectrum  of  roughness,  while 
7js(r)  is  the  short-wavelength  component,  also  generated  randomly.  The  short- 

wavelength  components  are  assumed  to  be  the  components  generating  scatter,  thereby 
responsible  for  reverberation,  predominantly  at  the  Bragg  wavenumber.  The  short 
wavelength  components  may  be  thought  of  as  the  easily  observeable  variability  of  the 
seafloor,  i.e.  loose  impediments  and  small  mounts.  The  longer  wavelengths  are  then  the 
larger  undulations.  Hence,  it  is  reasonable  to  think  of  a  ‘critical’  or  ‘limiting’ 
wavenumber  separating  these  two  components.  This  critical  or  limiting  wavenumber  is 

assumed  to  be  approximately,  kUm  ~^-k0d*,  where  6C  is  the  critical  angle  of 
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propagation,  and  smaller  roughness  wavenumbers  (i.e.  longer  wavelengths)  are  assumed 
to  dominate  the  forward  scatter  properties  of  the  acoustic  field. 

Performing  a  Taylor  series  expansion  of  the  Heaviside  step  function  we  have 
defined  in  Eq.  (2.4),  we  can  then  re-write  the  original  Helmholtz  equation,  Eq.  (2.1),  as 

V2p  +  k2(n2+ju)p  =  0  (2.6a) 

where 

n2  =nl(l-H(z-zb-r]l))+n2bH(z-zb-rjl)  <2-6b> 

and 

ju(r,  z)  =  (n2  -  n2  )rjsS(z  ~zb -?],).  ( Z6c ) 

This  shows  that  the  primary  part  of  the  propagation  is  due  to  the  long-wavelength 
components  and  the  perturbation  is  due  to  the  small  wavelength  components.  The  small- 
scale  features  previously  expressed  as  r/s(r)  are  considered  to  introduce  secondary  cross¬ 
terms  whose  contribution  will  be  small  and  may  thus  be  ignored.  From  hereon,  ignoring 
the  short-wavelength  components  will  introduce  an  approximation  to  the  theoretical 
development  and  variables  will  then  be  “hatted”  to  denote  the  use  of  only  long- 
wavelength  perturbations,  e.g.,  the  approximate  refractive  index,  n . 


We  now  express  the  wave  equation  in  cylindrical  coordinates  as 

lp  +  k2(n2+p)p  =  0. 


~ld_ 

d  ) 

i  a2  a2' 

r  dr 

1>J 

IV  a^2V2_ 

(2.7) 


Now  let  p  =  — prM ,  then 
Vr 


d2u  .  1  d2u  .  d2u  .  l2/„2 


•  + 


dr2  r2  d(p2  dz 


+^  +  k^  +  ju)u  =  0. 


(2.8) 


Next,  we  define  the  Q  operator  as 
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„  .2  i  id2  a2  ] 

Q-  n  +M+T2  ~2j72+TT  • 
A  kQ[r  dtp  dz 


In  general,  the  azimuthal  coupling  is  ignored  in  a  Nx2D  model  by  neglecting  the  — -r 

d<p 

term.  The  influence  of  azimuthal  coupling  may  be  expected  to  be  of  second  order,  so  the 
neglect  of  such  influences  here  is  consistent  with  the  treatment  of  small  perturbations. 


The  wave  equation  now  becomes 


i+k2Q2u=0. 


(2.10) 


Defining  u  in  terms  of  its  incoming  and  outgoing  components,  Tappert  (1977)  showed 
that  we  may  write  it  in  the  form 


=  — L(k+ +h_)  ,  and  ^~  =  ik0JQ(u+  -u_)  . 

V<2  dr 


(2.11) 


The  wave  equations  defining  the  evolutions  of  the  outgoing  and  incoming  fields  can  then 
be  shown  to  satisfy 


du+  ..  .  1  dQ 

-i-~ik0Qu++-—u_  ,  and 
dr  2  or 


d u .  1  dQ 


(2.12) 


This  shows  how  the  fields  form  a  set  of  coupled,  one-way  parabolic  wave  equations. 

Since  the  employment  of  an  operator,  Q,  is  not  algebraic  and  therefore 
inconvenient,  we  invoke  the  wide-angle  parabolic  approximation  (Thomson  and 
Chapman,  1983)  and  define  Qwape, 


Q  -  Qwap,  -  S+  ~L  ^  +[(«2  i]. 


(2.13) 


We  can  re-write  it  in  terms  of  a  Hamiltonian  differential  operator,  Hv 


Qwape  ~  1  (^  +  —  1  HWAPE(r, z,  -  2  )  » 


(2.14a) 
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where, 


H 


WAPE 


12  \ 


r,z. 


dz 2 


=  T  +  U 


is  written  in  terms  of  kinetic  and  potential  operators. 


T  =  l-  1+ 


1  d2 
kl  dz2 


and 


U  =  - 


(n2  +  juY2~  1 


Since  we  expect  pi  to  be  small, 


2  n  2  n 


(2.14b) 


(2.14c) 


(2.14d) 


(2.15) 


where  we  have  also  defined 


U(r,z)  =  -(n(r,z)~  1) 

as  the  unperturbed  potential  function.  Equation  (2.14a)  then  becomes 

d2 .  1  pi{r,z) 


Qwape~1  HwAPE(r,Z,dz2)+2n(r,z: ) 


(2.16) 


(2.17) 


since  HWAPE (r,  z,  ^  2 ) —  (T'  . 


i? 

'dz2 


Now,  ~  ^WAPE  ~  -Y~~  since  the  strongest  range  dependence  is  assumed  in 
dr  dr  2 h  dr 


pi.  Then,  the  evolution  of  the  outgoing  and  incoming  fields  from  Eq.  (2.12)  becomes 

(2.18a) 

(2.18b) 


du.  ..  .,  ,1  dpi  , 

ik0u+  —  ik0HWAPEu+  +  ^  ^  u_  ,  and 


.  a  1  $ju 

—  ~  -ik0u_  +  ik0HWAPEu_  +  u+ . 
dr  4 n  dr 

From  Eq.  (2.18a)  and  (2.18b),  the  two  fields  are  seen  to  couple  through  the  third  term  on 

the  right-hand-side  of  each  equation.  We  now  invoke  the  multiple-forward-single- 
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backscatter  approximation.  This  assumes  that  the  forward  propagation  can  only  scatter 
into  the  backward  (reverberant)  propagation  direction  once.  Since  energy  initially  only 
propagates  outward,  that  is,  only  u+  exists,  it  is  reasonable  that  this  will  contribute  to  u_. 

du 

This  is  seen  in  Eq.  (2.18b)  through  the  scattering  term  -r— .  With  the  single- 

dr 

backscatter  approximation,  we  then  say  that  the  secondary  scatter  influence  on  u+  by 
du 

— u_  in  Eq.  (2.18a)  is  negligible.  Thus,  the  forward  propagation  is  assumed  to  be 
or 

unaffected  by  the  small  scale  perturbations.  The  influence  of  fJ.  is  only  to  generate  source 
terms  for  the  backscattered  field. 


1.  Forward  Propagation 


We  can  now  obtain  u+  by  defining  u+  =  y/eik°r  so  that  Eq.  (2.18a)  (ignoring  the 
du 

—u_  term)  becomes 
dr 


dw 

+ ik0HWAPEy/  =  0 , 


(2.19) 


which  defines  the  evolution  of  the  forward  propagating  PE  field  function  if/.  The  solution 
to  this  equation  is  obtained  from  the  MMPE  model  with  the  long  wavelength  fluctuations 
built  into  n2 . 


2.  Backward  Propagation 

In  similar  fashion  to  the  forward  propagation  problem,  we  define  an  incoming 
field  function,  u_  =  <t>e~‘k°r ,  so  that  Eq.  (2.18b)  now  becomes 

=  -ik0< D + J-|^e‘Vyg<V  s  (2.20a) 

dr  4 h  dr 

noting  that  the  second  term  on  the  right-hand  side  is  the  forward  coupling  term  previously 
derived.  The  forward  coupling  term  /,  is  identified  as  the  inhomogeneous  term  in  our 
solution  for  <I>, 


•  10 


(2.20b) 


j  -  ^  jfyf  vrpik°<'r'*r) 

4nd/W 

Recall  that  y/  and  <!>  are  the  field  functions  of  the  two-way  propagation  with  the 
long  wavelength  components  built-in.  The  inhomogeneous  term,  I,  then  defines  the 
coupling  from  scattering  due  to  the  short  wavelength  perturbations.  In  the  next 
subsection,  we  will  introduce  the  point  source  and  then  cast  the  incoming  and  outgoing 
fields  into  the  form  of  the  sonar  equation  and  see  the  two-way  coupling  contributions 
likened  to  transmission  losses  but  with  a  wave  scattering  strength. 

3.  Forcing  Function 

Up  till  now,  we  have  only  dealt  with  the  homogeneous  wave  equation 

V2p  +  k2p  =  0,  (2.21a) 

with  no  hint  of  a  source  function.  If  a  point  source  is  now  introduced,  Eq.  (2.21a) 
becomes 

V2p  +  k2p  =  -4^P08(r  -rs),  (2-21b) 

where  P0  is  the  source  pressure  at  a  reference  distance  .  This  may  be  expressed  in 
terms  of  the  Green’s  function,  g(r,rs) ,  which  is  the  transfer  function  from  a  point  source 
to  a  point  receiver  defined  by  (Jensen,  et  al.,  1994) 

(y2+k2)g(r,rs)  =  -S(r-rs),  (2.22a) 

such  that  the  pressure  is 

p(r,  rs )  =  47tR0P0g(r,  rs  ) .  (2.22b) 

Extending  this  argument  to  an  extended  source,  we  can  apply  the  idea  of  superposition  of 
a  number  of  point  sources  to  say  that 

(V2  +  k2)p(r)  =  f(r)  =  -47t^w(rs)8(r  -  rs)d3rs ,  (2.23a) 

then 


11 


p(rs)  =  4^1  w(rs)g(r,rs)</3rs . 


(2.23b) 


We  can  also  define  the  Green’s  function  alternately  in  terms  of  the  PE  field 
function  ^by  letting 


4^g(r,Fs)  =  J%(r,rs)eik^. 
\  r 


(2.24) 


The  corresponding  inhomogeneous  PE  field  function  (Eq.  (2.19))  for  the  outward 
propagating  field  due  to  a  point  source  becomes 


^  +  iKHWAPE¥  =  ~4x^e*°rS(r  -  rs ) 


(2.25) 


since  \j/\  =  1  is  normalized  at  r  =  R0 .  Similarly,  the  inhomogeneous  form  for  the 


scattered  (or  backward  propagating)  field  is 


d<J> 

—  -ik0HWAPE®  =  -47& 


(2.26a) 


where  the  inhomogeneous  source  term  is  due  to  the  scattering  of  the  forward  propagating 
field,  defined  by 


3  = - —^~y/e,ko(r+r') , 

16 mdr 

which  acts  like  a  point  source  function  locally  at  the  scattering  point. 


(2.26b) 


We  now  need  to  define  some  geometrical  setup  to  aid  in  the  understanding  and 
development  of  the  effects  of  scatter.  Defining  the  location  of  the  scattering  point  as  rs 

and  the  volume  of  the  scattering  region  as  Vs ,  then  the  received  scattered  (backward 
propagating)  field  at  location  rR  will  be 


P-(rR)  =  4^o|3(^)J.  ^jre,kolrR'rslg(rR,rs)d3rs . 


(2.27) 
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Since  the  Green’s  function  is  the  transfer  function  between  two  given  points,  we  can  say 
that  g(rR,rs)  is  the  pressure  field  Green’s  function  from  each  point  within  the  scattering 

region  to  the  point  receiver  location  at  rs .  By  reciprocity,  we  say  that 

g(rR,rs)-g(rs,rR).  (2.2S&) 

Also 


then 


g(rs,rR) 


1 

4/zKo 


P+(fs,rR), 


(2.28b) 


P-(rR)  =  - 


1 

-Idju 

16^  l 

[hdr'\ 

rs~rR 

(rr,rDW^"%3n 


(2.29) 


1 


,,  -J^J7  P+rP+*A  > 

16^?0  ’nor 


where  p+T  is  the  forward  propagated  field  from  the  transmitter  to  the  scattering  point  and 
p+R  is  the  forward  propagated  field  from  the  receiver  to  the  scattering  point  which  by 
reciprocity  is  the  same  as  p_R . 


1  dll 

We  now  examine  the  integrand  term  ——7  .  Since  h2  is  defined  in  terms  of  the 

n  or 


step  function, 


n2  =nl{l-H(z-zb-nl))+n2bH(z-zb-T)i)  , 


(2.30) 


then  with  H(0)  =  — ,  the  value  of  n2  at  the  interface  is  just  the  average  value, 
2 


n2(z  =  zb+rjl)=]-(nl+nl). 


(2.31a) 


Solving  for  n  gives 
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(nl+n2) 


(2.31b) 


0  *  (nl  ~ nb ~zb~Vi), 


(2.32) 


where  we  have  ignored  the  slower  range  variations  in  zb ,  tj,,  nl  and  n2b  over  the 
scattering  region.  Substituting  the  expressions  in  Eq.  (2.31b)  and  (2.32)  into  Eq.  (2.29), 


we  obtain 


(2.33) 


V2  r  a?7,  2 
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where  p+Tb  and  /?+Ri)  are  the  values  evaluated  at  the  bottom  interface.  Note  that  the  range 
derivative  of  %  is  in  the  direction  from  the  scattering  patch  to  the  receiver. 

4.  Statistical  Treatment  of  the  Scattered  Field 

Of  interest  is  the  mean  squared  reverberation  from  different  roughness 
realizations,  which  from  Eq.  (2.33)  is  defined  as 


2 _ \nw~nb 


-  JJ(  • >  VM 


(2.34) 


V'-'  /  (16^)2  nl  +  n2b  JJ\|3r  dr'\ 

where  rs  and  rs'  indicate  two  different  boundary  realizations.  Consider  the  statistical 
characteristics  of  the  roughness.  We  note  that  the  small  scale  roughness  must  have  zero 
mean,  i.e.  (tjs)  =  0  and  has  an  amplitude  spectrum  of  S%  ZD(K,L )  for  a  given  realization, 
such  that 


Vs  (*,  y)  =  J  J  5%>2D  (*,  L)eiKxeiLydKdL , 


(2.35) 
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where  K  and  L  are  the  horizontal  wavenumbers  in  the  x-  and  ^-directions,  respectively. 
Here,  the  amplitude  spectrum,  S%2D(K,L) ,  refers  to  the  2-D  roughness  of  only  the 

small-scale  perturbations,  therefore  the  subscripts  ‘ rjs,2D  \  The  rms  roughness  is  then 


given  by 


rjs  |2 }  =  JJ  dKdl\\  dK'dL'  SVs,2D  (K,  L)S^2D  (K\  L')ei{K~K')xei{L-L'h 


Assuming  a  wide-sense-stationary  spectrum  in  space,  then 


(2.36) 


(|>7.f }  -  JJ dKdLWnaD  (K,L),  (2'37) 

where  WVs2D(K,L)-^S%2D(K,L)^  is  the  power  spectrum.  The  spatial  correlation 


function  for  the  roughness  can  then  be  defined  as 


(<7,  (*,  yH(x-x',y-  ft)  =  ||  dKdL  W,t2D(K,  De^e^ 


(2.38) 


Hence  for  a  given  realization, 


=  II k)WVtMe,K^%^dKdL , 


where 


(2.39a) 


kr=K  cos^  ,  (2.39b) 

and  </)  is  angle  in  the  direction  of  r  towards  the  receiver  relative  to  the  (x,y)  coordinate 
system.  The  following  diagram  shows  the  geometry  being  developed  here. 


Figure  2.1-  Scatter  Geometry 


Conversely,  the  alternate-domain  realization  may  be  expressed  as. 


IK 


dVs(x,y)  dt]s(x-x,y-y) 


dr 


dr 


)e'iK{x-xViL{y-y')d(x  -  x)d(y  -  y) 


(2.39c) 


=  *X.2  o(K,L). 


The  above  equations  may  be  considered  as  the  transform  pairs  for  achieving  a  bottom 
realization. 


Let  us  now  set  up  a  spatial  representation  of  the  notations  that  will  be  used  in  the 
ensuing  arguments.  Consider  two  scatter  points,  rs  and  r's  located  at  their  respective 

distances  from  the  transmitter,  T,  and  the  receiver,  R.  These  are  shown  in  the  following 
diagram. 


The  backward  propagation  or  scattered  field  in  Eq.  (2.34)  may  now  be  written  as 


drjs(x,y)  dr]s(x, y) 

1  1 

drR  drR 

4rTrR  A lrTrR 

,  ,  Wn  (rT )y/Rb  (rR (rT  )WRb  (rR  )  (2.40a) 


xeik’>(rT+r*~rhr!'))d2rsd2r' 


where 


(2.40b) 
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We  have  thus  far  assumed  that  the  fields  vary  slowly  in  range,  and  over  a  given 

scattering  patch  the  * —  amplitude  terms  also  vary  slowly.  Furthermore,  since  rjs  and 

y]rTrR 

decorrelate  rapidly  for  rT  £  rT  and  rR  £  r'R ,  then  we  may  assume  that  these  are  the 
dr 

same  in  the  amplitude  terms  and  the  locations  of  evaluation.  Hence,  the  field  functions 
become  independent  of  the  integral  over  the  scattering  patch  so  that 


(2.41) 


jy< 


dTjs(x,y)drjs(x,y') 


drR  drR 


eik«irT+rR-rl-rR) 


where  rT  =—  (rr  +  4),  rR=  —  (rR  +  r'R )  are  the  mean  values  of  the  distances  from  the 
2  2 

scattering  patch  to  the  receiver  and  transmitter. 


Introducing  difference  variables,  r"  =  r2  -  rT  and  rR  =  rR  -  rR ,  and  assuming  the 
roughness  is  isotropic,  WVs  2D(K,L )  =  WVs  2D(k) ,  we  can  write  the  above  integral  as 


drjs(x,y)  drjffx-x^y-/) 


drR  drR 


-**«+r*'d2  rsd2r's  =  J  KWVS'2D(k)d2rs  (2.42) 


where  k  =  &02 [(cos (j)T  +  cos <t>R)2  +  (sin <z>r  +  sin )2 J .  We  now  assume  that  over  the  small 
scattering  patch,  k2  and  WVs  2D(k )  vary  slowly  so  that 


<2'43) 

T  R 

and  A \  is  the  area  of  the  scattering  patch  on  the  bottom  interface.  The  new  geometry 
based  on  the  mean  distances  from  the  scattering  patch  is  shown  in  the  following  figure. 
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Figure  2.3  -  Scatter  Geometry  Based  on  Mean  Distances  from  Scatter  Patch 


To  simplify  the  analysis,  we  realign  the  diagram  to  have  the  receiver  line  coincide 

with  the  x-axis,  i.e.  (f>R  =  0.  Then  kr  =k  =■  k0cos<pT  and  k  =  2£0cos-^-.  Note  that  this 

wavenumber  indicates  that  the  spectrum  is  evaluated  at  the  Bragg  wavenumber.  The 
following  figure  shows  the  re-oriented  scatter  geometry. 


i 

Receiver,  rR 

V 

Transmitter, 

lid  Scatter  patch 

Figure  2.4  -  Re-Oriented  Geometry  for  <pR  =  0 


For  monostatic  scatter  then,  we  have 


(2.44) 
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where  r  =  rT=rR  and  |^|  =  |^n>l  =  W- 


Now  writing 


f 

nb=nw+An  =  nw 

\ 


\ 

) 


-  nw(\  +  8)  where  8  « 1 ,  then 


(2.45a) 


» (r2Sf 


=  n 


«  2nlS2(l-S)  ~  2n2S2  =  2(A nf. 


2(1+ S) 

_  ,  .  cn  c0  CoCc^-cJ  _  c0Ac 

Furthermore,  A n  =  nb-nw= - = - = - , 


so 


cb  c„ 


C  Cu 

wb 


C  Cu 
wb 


W-tM tot- 


C  Cu 
wb 


'Ac}2 

\c°  J 


,  since  cw~cfe~c0. 


Finally  for  the  monostatic  case,  the  backscatter  field  becomes 


v  c° ; 


(2.45b) 


(2.46) 


5.  Mean  Squared  Reverberation  Pressure  Level,  RPL 

Of  interest  is  the  mean  squared  reverberation  pressure  level,  RPL.  Reciprocity  of 
the  one-way  fields  has  been  invoked  in  order  to  write  the  solution  in  terms  of  the 
outgoing  field.  Using  the  solution  for  the  incoming  field  in  Eq.  (2.46),  we  can  now  cast 
this  into  the  form  of  a  SONAR  equation.  We  do  this  by  considering  the  contributions  of 
the  source  level,  directional  characteristics  of  the  source  and  receiver  and  the  area  of 
ensonification  being  captured  in  that  region  of  the  beams.  Hence  assuming  that  the  field 
functions  do  not  vary  rapidly  over  the  scattering  region  and  the  spectrum  of  the  roughness 
is  isotropic,  the  monostatic  reverberation  pressure  level  can  be  expressed  as 


( 

RPL  =  10  log 

V 


=  SL  +  DIt  +  DIR  + 10  log 


(2.47) 
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where  SL  is  the  source  level,  DIT/R  are  the  directivity  index  for  the  transmitter/receiver 
respectively,  and  A\  is  the  ensonified  area  of  the  bottom  interface  scattering.  The 

reverberation  loss  has  also  been  correspondingly  defined  as  a  general  function  of  either 
the  boundary  or  the  volume.  This  reverberation  loss  analysis  applied  to  the  boundary, 
from  the  analysis  in  the  previous  section,  is  defined  by 


where  Sb  is  the  full-wave  scattering  strength  due  to  the  small-scale  interface  roughness, 
and  WVs  2D(2k0)  is  the  2-D  spectrum  of  the  interface  roughness  evaluated  at  the  Bragg 
wavenumber. 
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One  should  note  the  following  points: 


•  If  the  statistics  of  the  interface  are  homogeneous  over  a  large  portion  of  the 
area  of  interest,  the  interface  reverberation  structure  is  controlled  by  two-way 
transmission  loss. 

•  Because  the  two-way  transmission  loss  is  dominated  by  forward  and 
backward  propagation  over  long-wavelength  features,  the  statistics  of  this 
signal  should  be  related  to  the  low  spatial  wavenumber-space  statistics. 

The  statistical  model  to  be  used  for  the  interface  roughness  will  be  treated  in  the 
next  chapter. 

B.  VOLUME  SCATTERING 

Having  trudged  through  the  analysis  of  the  bottom  interface  extensively,  the 
volume  scattering  development  will  take  on  the  same  general  form  with  the  exception 
that  the  index  of  refraction  varies  three-dimensionally  in  space.  We  will  specifically 
consider  for  now  only  the  theoretical  development  of  sound  speed  variations  and  the 
overall  effects  on  scattering  and  reverberation.  The  other  element  affecting  volume 
scattering  is  density  variations  in  the  volume,  which  will  not  be  treated  in  this  thesis  but 
will  be  left  for  follow-on  work. 

We  start  with  the  wave  equation  in  Eq.  (2.1)  and  the  index  of  refraction 
perturbations  in  the  spatial  scale  in  Eq.  (2.3), 

V2p  +  £02n2p  =  0,  (2.51) 

and 

n2  =n2wa-H(z-zb))  +  n2bH(z-zb),  (2-52) 

as  before,  but  now  the  bottom  (volume)  index  of  refraction  varies  three-dimensionally  in 
space,  nb(x,y,z ) .  As  before  in  the  statistical  treatment  of  the  interface  spatial  variation, 
sound  speed  variations  may  be  defined  by  (Yamamoto,  1995) 
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Cb  (x,  y,  z )  =  cb 0  (1  +  bz  +  8 (x,  y,  z)) , 


(2.53) 


where  cb0  is  the  mean  bottom  sound  speed  at  the  interface,  b  =  —  is  the  normalized 

Cb0 

gradient  of  the  bottom  sound  speed,  and  8  is  the  zero-mean  random  perturbation,  i.e. 
(£}  =  0  and 

(8(x, y, z),8(x, y', z))  =  JJjV*t30 (K, L,M)eiK(x~x)eiUy~y)eiM(z~z)dKdLdM  ,  (2-54) 

where  WS  3D(K,L,M)  represents  the  3-D  power  spectrum  of  the  fluctuations. 


Again  we  break  this  up  into  long  and  short  wavelength  components  by  having  the 


random  perturbations  split  as 

s=s,+s,, 

such  that  Eq.  (2.53)  now  becomes 

cb  (x,  y ,  z)  =  cb0  (1  +  bz  +  S, )  +  cb08s 
=  cb(x,y,z)  +  cb0Ss(x,y,z). 

This  can  be  incorporated  into  the  index  of  refraction  by  writing 


nb  = 


\C»J 


cb0  (1  +  bz  +  S,+  8S  )2 


»^a-2  Ss)  =  n2b+Mb, 


where  now 


(2.55) 


(2.56) 


(2.57a) 


jub=~2n;Ss  (2.57b) 

and  fil  is,  as  before,  the  description  of  the  environment  including  all  long  wavelenth 
component  perturbations. 

The  general  description  of  the  environment  of  Eq.  (2.52)  is  now 


n2  =  nxw(l-H(z-zb-  r)t ))  +  nbH(z  -zb-r}{). 


(2.58) 
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and  the  effects  of  interface  scatter  and  volume  scatter  may  be  considered  separately.  Note 
that  the  long-wavelength  contribution  to  the  interface  roughness  has  been  included  as  T]t 
into  the  perturbation,  which  in  turn  affects  the  effective  index  of  refraction.  The 
derivation  now  follows  exactly  as  before  arriving  at  the  final  result  for  the  bottom 
interface  in  Eq.  (2.33), 


P-~- 


f-i 


1  dfi„ 


P+tP+r  d  rs  > 


(2.59) 


16^Ro  £  h  dr 

where  the  volume  integration  occurs  over  the  sediment  volume  scattering  region,  i.e. 

Z>Zb+TJl. 


We  now  have 


^=*-2  n2 


dr 


Ml 

dr 


which  remains  a  volumetric  quantity,  so 


(2.60) 


r  „  dSs  ,3 

P-=T7-^)n-r: P+tP+r*  rs 


\67tRQ  l  dr 


(2.61) 


Taking  the  same  approach  as  the  interface  in  Eq.  (2.34),  we  see  that  the  mean  squared 
volume  reverberation  is  then 


(2.62) 


P+Tb P+Rb P+Tb P+Ro)d  rS^  rS  ' 

As  before,  we  write  the  pressure  fields  in  terms  of  the  field  function,  if/  (see  Eq.  (2.40a)), 
and  argue  that  they  do  not  change  significantly  in  range  within  the  scattering  region. 
They  do,  however,  vary  considerably  in  depth,  so 


=  To  1  \2  j~ ff  tyr (rr . ( rT , z)WR (rR,z)ys*R (rR , z)n(rs )n(r' )} 

'  '  (o^K0  j  rT  rR  zz 


(2.63a) 
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Recall  that  rT  =^-(rT  +  r')  and  rR  =  —(rR+r'R).  The  inner  term  is  evaluated  over 
difference  variables  (as  before  in  Eq.  (2.42))  r”  and  r"  to  give 


II 


^  (X,  y,  z)  (x -x",y- /, z  -  z)  \ 


yik°(r^4))d2rsd2n 


(2.63b) 


where  K",  L”  and  kr  are  as  before,  but 

W,sj2DB(£,L)^3D(tf,^  (2.63c) 

is  the  2D  spectrum  of  the  small  scale  perturbations  in  the  horizontal  dimension. 


If  we  now  assume  that  the  spectral  density  is  essentially  a  delta-correlated 
function  in  depth  and  horizontally  isotropic,  then  we  can  say  that 
wss,2d„  (K,L)  =  WSs2Dh  (k)  =  wss,3D(k,0)S(z  -  z) ,  and  neglect  variations  in  WSs2Dfi  over 
the  scattering  region.  Then,  assuming  a  monostatic  geometry, 

(W2)  =  -^j^-pr  ]  (n2(r,z)|^r(r,z)|2|^(r,z)|2)  (2-64) 

x^WSs  2DH(k)dZAA  , 

AA  being  the  horizontal  area  extent  of  the  scattering  volume,  much  akin  to  the  scattering 
patch  for  the  interface  reverberation  development. 

For  backscatter  then, 

(H2}=(— ^7-  J  (n2(r,z)|^r(r,z)|2|^(r,z)|2}  (2.65a) 

xk20WSs2Dii(2k0)dzAA. 

If  we  now  assume  that  the  spectrum  WSs  1Dh  does  not  change  significantly  with  depth, 
then 
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(2.65b) 


(k|2)4^;l  WSs'2DH(2k0)AA  ](n2(r,z)\Hr,z)\yz. 

1  V  J  z^i 

Finally,  we  will  assume  that  ^n2(r,z)|^(r,z)|4^«2^n(r,z)|^(r,z)|2^  ,  similar  to  our 
previous  treatment  of  fully  saturated  statistics. 


For  the  volume  reverberation  problem,  we  see  that  we  cannot  write  the  RPL 
simply  in  terms  of  transmission  loss  (TL)  but  instead  must  integrate  over  depth  the 

quantity  ^n(r,z)|^(r,z)|2^  at  each  range,  r.  We  shall  then  use  a  form  of  ergodicity  and 

write 


J ^n(r,z)|^(r,z)|2^  dz  ®  n(r,  z)|^(r,  zfdzj  ,  (2.65c) 

and  since  n  and  y/  are  based  only  on  the  long  wavelength  components,  we  will  simply 
say  that  we  can  use  (jn(r,z)|^(r,z)[2<iz)  .  Then  as  in  Eq.  (2.47),  we  may  write 


PPL  =  10  log 


„2 

Pref 

V  '  '  J 


=  SL  +  DIt  +DIr  +10  log 


-RLV  , 


(2.66a) 


where 


RLV=- 20  log 


■p  °°\h(r,z)\v(r,z)fdz 


z>zb+v, 


-10  log 


k2 

ZQ  w 

2  Ss,2D„ 


32x 


(2*n) 


(2.66b) 


In  the  above  expression,  2Z)(/ (2fc0)  is  the  2-D  horizontal  spectrum  of  the 

volume  fluctuations,  which  is  assumed  isotropic  and  independent  of  depth  (as  in  the 
interface  reverberation  problem)  evaluated  at  the  Bragg  wavenumber,  and  the  depth 
integral  extends  from  the  bottom  interface  downwards  to  negative  infinity.  Note  that 
volume  reverberation  loss  cannot  be  expressed  simply  in  terms  of  the  two-way 
transmission  loss  as  before. 
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One  should  also  note  the  following  few  points: 

•  The  reverberation  loss  is  not  independently  influenced  by  the  forward 
propagating  fields,  but  is  a  magnitude  scaling  or  weighting  function  of  the 
local  index  of  refraction  (positive  definite). 

•  If  the  statistics  of  the  volume  are  homogeneous  over  a  large  portion  of  the  area 
of  interest,  the  volume  reverberation  structure  is  controlled  by  this  non-linear, 
two-way  transmission  loss  between  sound  speed  structure  and  acoustic  field. 

•  Because  this  weighted  two-way  transmission  loss  is  dominated  by  forward 
and  reciprocal  propagation  over  long-wavelength  features,  the  statistics  of  this 
signal  should  be  related  to  the  low  spatial  wavenumber-space  statistics,  as  in 
the  case  of  interface  reverberation.  We  will  however  expect  to  see  smoother 
changes  (not  abrupt)  as  compared  to  the  interface  reverberation  statistics.  This 
is  because  the  statistics  apply  non-linearly  here  providing  only  the  mean- 
square  positive-definite  weighted  average  of  the  acoustic  interactions  with  the 
environment. 

The  statistical  model  to  be  used  for  the  volume  sound  speed  variations  will  be 
treated  in  Chapter  HI. 

C.  TIME-DOMAIN  ANALYSIS 

1.  Time-Domain  Analysis  of  the  Interface 

Up  until  now,  the  focus  has  been  on  CW  analysis.  The  previous  analyses  for  both 
the  interface  and  volume  have  suggested  that  the  structure  of  the  reverberation  is  due 
mainly  to  the  effects  of  propagation.  To  predict  the  effects  of  a  pulsed  signal  instead  of 
the  CW  signal,  the  MMPE/MMPEREV  model  is  run  over  a  spectrum  of  frequencies. 
Subsequently,  Fourier  analysis  will  give  the  time-domain  structure  of  the  pulse 
propagation  in  time. 
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In  shallow  water,  propagation  of  a  signal  from  a  point  source  will  experience 
multipath  effects.  In  the  time-domain,  a  scatter  patch  receiving  multipath  arrivals  will 
scatter  the  received  signal  at  different  arrival  times.  The  scattering  patch,  by  reciprocity, 
scatters  the  entire  received  signal  back  through  the  same  multipaths  (monostatic 
geometry  assumed)  that  each  signal  has  traversed,  leading  to  the  receiver  receiving 
multiple-multipath  signals  at  different  times.  Thus,  the  scattering  patch  is  a  source, 
producing  time-delayed  replicas  (with  proper  amplitude  weightings)  back  to  the  receiver 
through  the  multipath  structure  of  the  water  column.  The  geometry  of  the  development 
may  be  represented  by  the  following  diagram. 


As  the  source  transmits  a  signal  out  to  the  scattering  patch,  the  scattering  patch 
receives  signals  along  all  the  multipaths  (in  the  figure  above,  just  three  multipaths).  The 
scattering  patch  scatters  each  signal  and,  by  reciprocity,  these  scattered  signals  follow  the 
same  (three)  multipaths  back  to  the  source  location.  Recall  that  the  Green’s  function  is 
the  transfer  function  between  the  source  and  the  receiver  and  it  does  not  matter  which 
direction  the  propagating  wave  is  travelling.  Hence,  from  a  signal  point  of  view,  these 
may  be  represented  as  depicted  in  Figure  2.6.  The  total  reverberation  signal  from  this 
scattering  patch  would  then  be  the  coherent  sum  of  the  three  lower  panels. 


27 


>ignal  from  Path  1 


Path  2 


In  general,  if  the  receiver  and  source  are  not  co-located,  it  may  be  useful  to 
consider  the  arrival  times  of  the  same  signal  from  the  source  and  receiver  to  the  same 
scatter  point/patch.  By  reciprocity,  we  may  say  that  the  receiver  has  “propagated”  the 
signal  towards  the  scatter  patch.  The  geometry  has  been  adapted  from  Smith  and 
Cushman  (1997)  as  shown,  in  Figure  2.7. 


Figure  2.7  -  Two-way  return  from  a  Scattering  Patch 
These  two  propagations  may  be  considered  independent  and  in  order  to  compute  the  total 
reverberant  field  at  the  receiver,  we  may  convolve  the  two  field  functions  as  explained  by 
Smith  and  Cushman  (1997).  With  an  understanding  of  the  pulsed  propagation  and 
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multipaths,  the  one-way  time-domain  pressure  due  to  a  pulsed  signal  at  the  scattering 
patch  is 

P+T/Rb(rm’ty  =  J 

where  the  forward  propagated  pressure  field  from  the  transmitter/receiver  at  the  interface, 
p+T/Rb ,  is  specified  at  range  rm,  which  is  the  range  at  the  mh  range  step.  The  field- 

function  evaluated  at  the  interface  at  range  rm  and  at  frequency  /  may  be  extracted  from 
the  forward  propagation  of  the  broadband  field  computed  by  MMPE.  Notice  that  in  order 
to  arrive  at  the  required  pressure  in  the  time  domain,  we  have  to  integrate  over  all  the 
frequencies,  and  the  above  expression  is  really  the  Fourier  transform  of  the  field  in  the 
frequency  domain. 

Extending  the  argument  further,  the  travel  time  of  the  two-way  pressure  field  at 
the  receiver  is  now  the  convolution  (Smith  and  Cushman,  1997)  of  two,  one-way  fields  in 
the  time-domain  propagated  forward.  This  is 

P 2-way,, b  C Ym  ’ 7)  =  J  P+Tb  ( rm '  *)P+Rb  ’ 7  _  0^  *  (2.68a) 

where  we  recall  from  Eq.  (2.33)  that  p+Tb  and  p+Rb  are  the  forward  propagated  fields 

from  the  transmitter  and  receiver  to  the  scattering  point  evaluated  at  the  bottom  interface, 
respectively.  Note  that  the  forward  propagated  field  from  the  receiver  to  the  scattering 
point  is  by  reciprocity  the  same  as  p_Rb .  We  also  note  that  the  transmitter  and  receiver 

need  not  be  co-located  in  depth. 

Applying  the  above  approach,  we  see  that  the  time-convolution  of  the  two  field 
functions  is  also  the  frequency  multiplication  of  these  functions  in  the  frequency  domain. 
Hence,  when  the  forward  propagating  fields  are  already  available  in  the  frequency 
domain,  we  may  then  compute  the  two-way  field  in  the  frequency  domain  from  the 
interface  as 

P  2-way, b  (*m  >  f)  =  P+Tb  (rm  >  f)P+Rb  *  /)  *  (2.68b) 


4Z 


¥T/Rb(rm’f)e 


.27? 

: — —r 
Co 


-n^df. 


(2.67) 
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where 


P+n ( rm . /)  =  ~j=¥+n ( rm , />' V ,  and 
•Jr„ 


(2.68c) 


P+Rb(rm’f )  / ¥+Rb(rm’f)e  » 

vr- 

at  a  particular  bottom  depth,  zb{rm) ,  at  range  rm,  and  at  each  frequency.  By  applying  Eq. 


(2.67),  we  can  then  arrive  at  the  two-way  travel  time  structure  of  the  reverberation  loss 
( RLb )  due  to  a  single  bottom  patch,  i.e. 


P.(rm,  1)  =  Aj  <2'68d> 

where  the  constant  A  has  been  included  to  incorporate  all  the  other  factors  needed  to 
define  reverberation  loss. 


We  then  continue  our  propagation  through  the  entire  water  space  of  interest  and 
arrive  at  the  reverberant  field  at  each  range  step.  To  see  the  total  field  at  the  receiver,  we 
coherently  sum  up  all  the  pressure  values  from  different  range  segments,  rm,  by  matching 
up  the  discrete  arrival  times,  tn, 

P.fr.)-tp.(r..O  <269) 

m- 1 

where  p_  is  the  total  pressure  at  the  receiver  received  at  time  t„.  Note  that  this  matching 
requires  a  coordinated  scaling  between  range  step  size  and  time  step  size,  as  will  be 
explained  in  Chapter  IV. 


2.  Time-Domain  Analysis  of  the  Volume 

The  volume  analysis  is  a  2-D  problem  in  that  the  reverberation  is  not  directly 
proportional  to  the  two-way  transmission  loss,  but  requires  proper  weights  through  the 
entire  volume  when  computing  the  reverberant  field  at  the  receiver.  The  development 
will  follow  closely  that  of  the  previous  sub-section  on  the  interface. 

Firstly,  the  reverberant  field  due  to  each  depth/range  point  is  computed  using  the 

frequency-multiplication  of  the  source-patch  and  patch-receiver  (by  reciprocity,  the 
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forward  propagation  used)  propagating  fields,  and  then  applying  a  Fourier  transform  to 
the  time-domain.  This  is 

Pi-way  (rm,Z,f)  =  n{rm ,  z)p+T  (rm ,  z,  f)p+R  (rm  ,  z,  /) ,  (2.70a) 

where  now  the  two-way  reverberation  signal  is  computed  for  every  grid  point  of  interest 
(z  >  zb  always)  at  a  particular  frequency,  /,  and  n(rm,z)  is  the  local  index  of  refraction 
at  the  grid  point.  The  justification  for  multiplying  by  the  local  index  of  refraction  is  to 
provide  the  same  weighting  used  in  the  CW  analysis.  Fourier  transforming  gives 

Pi-vay  (rm  >Z’t)  =  j p2_way  (rm ,  z,  f)e~i2^df .  (2J0b) 

The  2-D  problem  may  now  be  reduced  to  a  1-D  problem  by  performing  a  vertical 
integration  to  arrive  at  a  single  reverberation  time-series,  i.e. 

CO 

P-(rm,t)  =  B  jp2_way(rm,t)dz,  (2.71) 

J>2t 

which  is  the  coherent  sum  over  all  depths  below  the  interface  at  range  step  m.  Note  that 
the  constant  B  has  been  introduced  to  account  for  all  remaining  terms  in  the  reverberation 
loss,  RLv.  We  have  now  ‘collapsed’  the  volume  in  2-D  to  1-D,  giving  a  single  set  of  time- 
series  which  can  then  be  matched  and  summed, 

(2-72) 

m=\ 

giving  a  single  time-series  in  p_(t„)  which  is  the  two-way  time-domain  field  defining  the 
volume  reverberation  loss  at  the  receiver  at  time  t„  due  to  the  entire  volume  in  the  space 
of  interest. 

In  the  post-processing  analysis  to  be  discussed  in  Chapter  IV,  we  will  implement 
the  above  discussions  for  both  the  interface  and  volume  to  arrive  at  the  time-domain 
analysis. 
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III.  GENERATING  REALIZATIONS  -  IMPLEMENTATION  OF 
THE  REVERBERATION  PROBLEM  IN  MMPE 


Having  developed  the  theoretical  basis  for  modeling  the  interface  and  bottom 
variabilities,  we  now  have  to  define  statistical  methods  or  models  to  generate  the 
interface  and  sound  speed  variations.  We  first  developed  the  theoretical  aspects  of  the 
models  used  for  generating  both  the  interface  and  volume  perturbations,  then 
implemented  the  algorithms  in  MATLAB  to  enable  checking  of  the  theoretical 
developments  before  integrating  them  as  functions  into  the  MMPE  model  (Smith,  2000) 
written  in  FORTRAN  95.  Much  of  the  theoretical  development  for  generating 
perturbations  to  the  interface  is  based  on  the  work  of  Goff  and  Jordon  (1988)  while  the 
model  for  the  volume  is  based  on  Yamamoto’s  (1995)  work. 

We  should  note  that  generating  good  realizations  would  form  an  important  first 
step  in  the  effort  to  model  shallow  water  reverberation  and  their  effects  on  sound 
propagation.  Of  particular  importance  are  the  “tweaking”  that  must  be  afforded  in  root- 
mean-square  values  when  generating  these  realizations  so  that  we  are  able  to  study  the 
effects  of  different  bottom  undulations  and  sediment  types. 

A.  SPECTRAL  DESCRIPTIONS  OF  BOTTOM  VARIABILITY 

1.  Interface  Roughness  (Spatial  Perturbation) 

From  Chapter  II,  we  have  seen  that  the  interface  effects  can  be  separated  into  both 
long-  and  short-wavelength  fluctuations.  We  have  principally  neglected  the  short 
wavelength  components  in  the  forward  propagation  problem.  Since  the  interface 
roughness  is  essentially  2-D,  we  can  assume  a  2-D  interface  spectrum  of  the  form 

- — 5/  .  *,=VF+7F,  (3-D 

a+O,2) /2 
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where  Lcorr  is  a  correlation  length  scale,  is  a  spectral  exponent,  and  kr  is  the 


horizontal  spatial  wavenumber  vector.  Note  that  for  kr  » 


wV5'2D(kr»Lcorr)~k;‘ 


The  normalization  factor  ju  is  defined  in  terms  of  the  root-mean-square  (rms)  roughness 
a2  by  requiring 


27rjWr]s2D(kr)krdkr=(J2,  (3.3) 

o 

which  leads  to 

°Aa> 

n2  J 

Note  that  kr  is  the  horizontal  component  of  the  wavenumber,  and  the  2-D  spectrum, 
WVs  2D ,  is  assumed  isotropic  (independent  of  direction). 


For  this  description,  the  exponent  is  constrained  by 

2<j3v<4.  (3.4b) 

Consistent  with  previous  reverberation  work  by  Goff  &  Jordan  (1988)  on  a  rough 
interface,  we  shall  assume  a  fractal  character  with  /?  =  3.5 . 

For  the  scattering  amplitude,  we  may  simply  evaluate  WVs  2D{kr)  at  the  Bragg 
wavenumber,  kr  =  2k0.  However,  for  the  long- wavelength  interface  roughness,  we  need 
the  1-D  spectrum  along  the  x-axis,  WVs  W(K) .  We  do  this  by  taking  the  1-D  transform  of 
W%  2d(K,  L)  along  the  slice  y  =  0 , 

WVS'1D(K)=  ]w%'2D(K,L)dL.  (3.5a) 
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In  cylindrical  coordinates,  we  have 


where 


Tk rt^w2(kr)dk. 


A.  i 


=  ytT2L_(l  +  C^2)  2  2  , 


r  = 


(B„  \ 

^  - 1 

M  i) 

2  11 

1 11 

[ 2  hi 

m 


v2; 


(3.5b) 


(3.5c) 


We  now  have  to  transform  a  random  realization  of  the  corresponding  1-D 
amplitude  spectrum  which  has  been  scaled  by  a  random  amplitude  and  phase.  In  other 
words,  we  define  the  realization  as 


m=  fs^K^dK  (36a) 

— oo 

where 

S,(K)  =  [wVS'W(K)Y2A(K)ei6(K) ,  (3.6b) 

and  A  and  0  are  random  numbers  for  all  values  of  K.  Since  the  complex  amplitude  of 
each  wavenumber  component,  Ae,e ,  should  exhibit  a  normal  distribution,  the  random 
phase  of  each  component  can  be  obtained  from 

6  -  2mx ,  (3.7a) 

where  r\  is  a  uniformly  distributed  random  variable  in  the  interval  [0,1].  The  magnitude 
A,  however,  exhibits  a  Rayleigh  distribution.  Because  we  are  in  fact  generating  a 
realization  of  the  power  spectrum,  we  must  consider  the  magnitude-squared,  which  has  a 
negative  exponential  distribution.  The  random  amplitude  of  each  wavenumber 

component  is  then  obtained  by 
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A  =  V-  ln(r2),  (3.7b) 

where  r2  is  another  independent,  uniformly  distributed  random  variable  in  the  interval 

[0,1]. 


In  practice,  we  could  simply  use 


A,  i 


W,„w(.K)  =  (l+L2„„K2)i ** 


(3.8a) 


T](X) 


and  rescale  the  result  by  its  rms  value,  i.e.  — - - tt  ,  or  we  could  simply  say 


w. 


Vs 


m~(i 


+  L; 


2 

'corr 


£2)"^  2 


(3.8b) 


where  Lcorr  and  had  to  be  specified.  This  was  the  generic  spectral  model  used  in 

generating  the  realizations  first  in  MATLAB  and  then  for  implementation  in  the  MMPE 
model. 


2.  Volume  Sound  Speed  Fluctuations 

Unlike  spatial  variations  of  the  interface,  volumetric  variations  are  generally  3-D 
in  nature.  We  will  now  treat  the  generation  of  volumetric  variations  through  spectral 
realizations  as  in  the  interface.  The  sediment  volume  sound  speed  perturbations  may  be 
modeled  by  a  3-D  volume  spectra  (Yamamoto,  1995)  given  by 

W^d(K,L,M)  =  ^^-(a2(k2  +  L2)+M2)'^'1  (3-9) 

where  B  is  the  spectral  strength  constant,  J3S  is  the  spectral  exponent,  and  A  =  —  =  —  is 

a\  a2 

the  horizontal-to-vertical  aspect  ratio  describing  the  anisotropy  of  fluctuations  in  the 
sediment.  We  shall  assume  the  major  axis  corresponding  to  M  is  aligned  with  the  vertical, 
thus  we  neglect  any  “dip”  angle  in  the  structure.  M  is  thus  the  vertical  wavenumber. 
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From  Yamamoto’s  findings,  we  shall  assume  the  following  values: 


A~ 2, 

B  ~  5xl0-4 ,  (310) 

A  ~  5 . 

This  should  produce  a  relative  velocity  variance  of  ~10'4.  Also  note  in  this  form  that  we 
really  assume  K  =  \kx  where  fy,  =lm  is  simply  a  reference  length  so  that  K,  L,  M  are 

unitless  and  units  of  W^  3D  do  not  depend  on  ($. 


Recalling  Eq.  (2.66b),  in  order  to  evaluate  the  scattering  due  to  the  volume 
perturbations,  we  need  a  representation  for  the  volume  spectra  in  the  2-D  horizontal 
plane.  Hence  for  our  reverberation  computations,  we  have 


OO 

CST.L)  =  jwtsM(K,L,M)dM 


(3.11) 


Note  that  we  are  effectively  evaluating  the  3-D  spectra  at  a  single  depth,  since  it  is  delta- 
function  correlated.  Substituting  Eq.  (3.9)  into  (3.11),  we  obtain, 


OO 

WSs'2Dh(K,L)=  jws$'3D(K,L,M)dM 


(3.12) 


= MTrfr  (k2 + l2  )+  m 2  ]" 

7T  » 


7F 


rfft+il 

l 2  2 

2;r 


2  2  — X - 


+  1 


For  Ps-2,  Eq.  (3.12)  can  also  be  shown  to  reduce  simply  to 


wSM(K,L)=— |[A2(js:2+g)+M2]'!dM 

^  n 


=^[a2(^+l2)]-!. 


(3.13) 
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For  the  values  of  B  and  T  in  Eq.  (3.10)  chosen,  this  reduces  to 

»'Js.2c,(JC,L)  =  ff(jS:J+tJ)rl  (3.14a) 

where 

a  =  —  =  5xl0"5.  (3.14b) 

2A 

Note  that  the  above  provides  the  2-D  horizontal  spectrum  which  is  needed  to 
evaluate  the  reverberation.  However,  to  compute  the  forward  propagation,  we  need  the  2- 
D  vertical  spectrum  in  the  (r,z)  plane.  In  other  words,  we  will  now  perform  the 
integration  cross-range  at  y  =  0 ,  giving 


(3.15) 


(3.16a) 


(3.16b) 


To  generate  2-D  vertical  volume  fluctuation  realizations,  we  do  essentially  the 
same  thing  as  in  the  case  of  the  interface.  Since  we  are  looking  at  sound  speed 
perturbations,  we  define  a  realization  as 


c0S(x,  z)  =  JJ  SSs2Dv  (K,  M  )eiKxemd.KdM 


(3.17a) 
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where 


Ss„^K,M)A^s,m(K,M)Y1A(K,M)emKM).  (3.17b) 

Notice  that  since  we  have  treated  sound  speed  perturbations  in  the  volume  in  the  vertical, 
we  really  are  generating  a  series  of  vertical  realizations  at  each  range  step. 

In  similar  fashion  to  the  interface,  we  define  the  2-D  random  phase  variations  as 
d{K,M)  =  2jtrl{K,M)  ,  (3.18a) 

where  rx(K,M)  is  now  a  matrix  of  uniformly  distributed  random  numbers  in  [0,1],  and 
the  amplitude  variations  are 

A(K,M)  =  J^jjUA))  .  (3>18b) 

where  r2(K,M)  is  another  matrix  of  uniformly  distributed  random  numbers  in  [0,1].  In 
practice,  we  use 

WSim(,K,M)~{/?K2 ,  (3'19) 

where  A  and  f}s  are  to  be  specified  in  generating  the  realization.  The  final  result  for  the 
volume  perturbation  is  then  rescaled  by  the  desired  rms  value. 


3.  Spectral  Filter 

We  must  note  that  both  the  interface  and  volume  realizations  will  need  to  be 
filtered  at  low  and  high  wavenumbers  in  order  to  reduce  edge  effects,  which  will 
introduce  additional  unwanted  spectral  components.  The  filter  is  also  used  to  extract  the 
pertinent  portions  of  both  the  interface  and  volume  spectra  for  creation  of  the  long- 
wavelength  realization  affecting  the  forward  propagation. 


At  the  low  end,  we  can  assume  scales  greater  than  500m  are  known.  As  such. 


min 


2  n_ 
500 


=  0.1 


(3.20a) 
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At  the  high  end,  we  should  compute  up  to  about  the  acoustic  wavelength,  so 


which  will  depend  on  our  frequency  of  calculation.  In  a  broadband  computation,  c0  will 
be  our  center  frequency. 

B.  IMPLEMENTATION  OF  THE  INTERFACE  AND  VOLUME 
SPECTRAL  DESCRIPTIONS  IN  MMPE 

As  the  MMPE  is  well  documented  (Smith,  2000),  a  quick  run-through  will  be 
afforded  in  order  to  understand  the  integration  of  the  spectral  perturbation  model  into  the 
MMPE.  However,  before  any  implementation  was  possible,  MATLAB  was  used  to  test 
the  model  and  its  algorithm. 

1.  Description  of  MMPE 

The  N x2D  MMPE  model  (Smith,  2000)  is  based  on  the  split-step  Fourier 
algorithm  developed  in  the  UMPE  model  (Smith  and  Tappert,  1994)  written  in 
FORTRAN  ‘95.  MMPE  generates  the  environmental  propagator,  Uop,  in  a  subroutine 
called  “ENVPROP”  using  the  bottom  bathymetric  and  sound  speed  data.  These  data 
points  are  interpolated  to  fit  the  depth  points  being  computed  at  each  range  step. 
ENVPROP  is  subsequently  called  at  each  range  to  generate  the  propagator  variations  due 
to  a  range-dependent  environment  specified  in  the  input  files.  In  the  case  where  there  are 
perturbations,  the  program  is  ‘forced’  to  think  that  it  is  in  a  range-dependent 
environment. 

Several  input  files  are  used  to  specify  pertinent  input  parameters,  most  important 
of  all  being  the  file  named  “pefiles.inp.”  In  there,  several  input  filenames  are  specified, 
one  each  for  the  source  parameters,  sound  speed  profile(s),  bathymetric  data  set(s), 
bottom  properties  (like  shear  and  compressional  sound  speeds,  etc.),  a  sub-layer 
bathymetry,  and  the  sub-layer  properties.  A  total  of  seven  input  files  are  used  for  running 
MMPE  simulations. 
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2.  FORTRAN  &  MATLAB  Implementation  of  Perturbation 
Models 

Before  actual  implementation  of  Eqs.  (3.8b)  and  (3.19)  into  the  MMPE’s 
FORTRAN  code,  MATLAB  was  used  to  code  up  and  test  the  algorithm  and  effects  of 
generating  different  roughness  and  perturbation  rms  values.  These  were  proven 
successful,  with  the  expected  spatial  and  sound  speed  perturbations  modeled  effectively. 
A  “tweaking”  factor  for  each  type  of  perturbation  was  used,  namely  the  rms  roughness 
and  the  rms  sound  speed  variations,  to  allow  observation  of  these  perturbations. 

Subsequently,  these  algorithms  were  translated  and  implemented  into  the  MMPE, 
now  renamed  “MMPEREV.”  In  MMPEREV,  two  additional  subroutines  are  called 
within  the  subroutine  “ENVPROP1”  to  introduce  perturbations.  “ZBGEN”  introduces  the 
spatial  perturbation  to  the  bottom  bathymetry  at  each  range  step  while  “DCGEN” 
introduces  the  sound  speed  perturbation  to  the  volume  at  each  depth  grid  point  and  range 
step.  ENVPROP1  is  then  made  to  think  that  a  range-dependent  case  is  present  simply 
because  the  perturbations  cause  irregularities  in  depth  and  sound  speed,  causing  the 
program  to  call  ENVPROP  at  each  range  step. 

A  random  number  generator  was  used  in  FORTRAN  to  provide  the  uniformly 
distributed  random  numbers  in  the  range  [0,1].  Each  time  the  random  number  is  required, 
a  series  is  generated.  However,  each  time  the  compiled  MMPEREV  program  is  started, 
the  same  seed  is  used  by  FORTRAN  so  that  the  same  set  of  random  numbers  are 
invoked.  Thus,  the  same  sequence  of  random  numbers  is  used  for  each  run,  allowing  for 
controlled  analysis  of  the  results.  This  is  described  later.  Later  implementations  to 
MMPEREV  will  allow  for  the  generation  of  a  random  seed  in  order  to  effectively  model 
truly  random  fluctuations  to  the  intended  data  set.  Some  results  of  the  controlled  analysis 
are  seen  in  the  next  sub-section. 

3.  Results  of  Varying  RMS  Values  to  Perturbations 

We  now  present  some  results  from  the  random  fluctuations  to  bottom  bathymetry 
in  range  by  varying  the  root-mean-square  values  of  the  interface  spatial  perturbations. 
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Depth,  m  Depth, 


Interface  rms  roughness:  1m 


0  0.5  1  1.5  2  2.5  3  3.5  4  4.5  5 

Range,  km 

interface  rms  roughness:  5m 


0  0.5  1  1.5  2  2.5  3  3.5  4  4.5  5 

Range, km 

Figure  3.1-  Plot  of  Interface  Roughness  of  Different  rms  Values 


Interface  rms  roughness:  10m 


From  Figure  3.1  and  Figure  3.2,  increased  rms  roughness  is  seen  as  the  root- 
mean-square  values  are  increased.  The  longer  wavelength  perturbations  also  become 
more  distinct.  A  cursory  survey  of  shallow  water  areas,  especially  in  the  East  China  Sea 
where  the  ASIAEX  experiment  will  be  conducted,  reveals  that  the  long  wavelength 
components  are  relatively  small.  Hence,  an  rms  value  of  lm  was  selected. 

Effects  of  tweaking  volumetric  root-mean-squared  perturbations  to  the  sound 
speed  follows.  The  mean  bottom  sediment  sound  speed  selected  was  1700  m/s,  which  is 
normal  for  a  sand/mud  layer. 

Sound  Speed  Profile  Data  Over  Depth  &  Range,  (m/s),  Volume  Perturbation:  5m/s 

1710 

1700 

1690 

1680 

1670 

1660 

1650 

1640 

0  0.5  1  1.5  2  2.5  3  3.5  4  4.5 

Range  (km) 

Figure  3.3  -  SSP  Data  with  Volume  Sound  Speed  Perturbation  of  5  m/s 
and  Interface  Perturbation  of  1  m 


43 


Depth  (m)  hcj  Depth  (m) 


Sound  Speed  Profile  Data  Over  Depth  &  Range,  (m/s),  Volume  Perturbation:  15m/s 
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1700 


1680 


1660 


1640 


Range  (km) 

igure  3.4  -  SSP  Data  with  Volume  Sound  Speed  Perturbation  of  15  m/s 
and  Interface  Perturbation  of  1  m 


Sound  Speed  Profile  Data  Over  Depth  &  Range,  (m/s),  Volume  Perturbations:  25m/s 


Range  (km) 


Figure  3.5  -  SSP  Data  with  Volume  Sound  Speed  Perturbation  of  25  m/s 
and  Interface  Perturbation  of  1  m 
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Apparent  from  Figure  3.3  through  Figure  3.5  is  the  increasing  contrast  in  the  sound  speed 
profile  for  an  increasing  sound  speed  rms  perturbation  when  comparing  them 
qualitatively.  While  it  is  not  likely  that  the  sediment  layer  is  homogeneous,  as  in  Figure 
3.3  (rms  perturbation  being  5  m/s),  it  is  not  unimaginable  that  the  sound  speed  profile  can 
vary  considerably  due  to  inhomogeneities  and  mixture  of  different  sediment  types.  A 
median  sound  speed  perturbation  of  15  m/s  (Figure  3.4)  was  finally  used  for  the  follow- 
on  analysis.  Notice  also  that  the  interface  rms  roughness  used  was  1  m. 
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IV.  POST-PROCESSING  IMPLEMENTATIONS  IN  MATLAB  AND 

SIMULATION  RESULTS 


The  second  and  final  step  in  our  modeling  efforts  was  the  development  of  post¬ 
processing  routines  to  extract  the  relevant  data  from  two,  one-way  propagation  data  sets 
generated  by  MMPEREV  (see  Chapters  II  and  El)  and  then  compute  reverberation  and 
perform  signal  analyses  on  the  derived  data.  We  will  also  see  that  the  time-domain 
convolution  of  the  two  one-way  field  functions  will  give  the  time-domain  reverberation 
structure  received  for  both  monostatic  and  bistatic  setups.  Routines  were  written  to 
perform  signal  processing  on  these  reverberation  data  sets.  Specifically,  this  chapter  will 
expand  on  the  work  done  to  derive  the  following: 

•  reverberation  loss  for  broadband  signals  (time-series)  and  CW  signals  (range- 
series); 

•  vertical  correlation  and  peak  correlation  computation;  and 

•  spectral  analyses. 

We  saw  from  Chapter  II,  Sections  B  and  C  that  a  scaling  or  normalization  due  to 
the  wave  scattering  strength  is  needed  in  order  to  arrive  at  the  reverberation  loss. 
However,  we  will  not  attempt  to  introduce  any  scaling  here,  as  these  normalization 
factors  are  assumed  to  be  constants  (see  Eqs  (2.50)  and  (2.66b)  for  CW  analyses,  and 
Eqs.  (2.68b),  (2.68c)  and  (2.71)  for  time-domain  analyses). 

Discussion  of  the  reverberation  computation  cannot  be  made  without  first 
defining  the  geometrical  and  environmental  setup  of  the  modeling  effort.  Of  interest  was 
the  particular  configuration  and  environmental  parameters  likely  to  be  used  in  ASIAEX. 
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A.  MODELING  ENVIRONMENT  AND  GEOMETRY 


Since  the  varying  spatial  properties  of  the  broadband  signal  are  of  interest,  an 
array  was  typically  chosen  in  collecting  field  data.  Each  element  of  the  vertical  line  array 
(VLA)  was  assumed  to  be  a  point  source  and  the  array  was  located  vertically  in  the  water 
column.  The  MMPE/MMPEREV  model  provided  as  an  output  the  field  function  for  the 
entire  water  column  and  bottom  volume  up  to  a  given  range.  The  fluctuations  in  the 
bottom  bathymetry  and  volume  sound  speed  were  randomly  generated  by  MMPEREV  as 
explained  in  Chapter  HI. 

1.  Multi-static  Reverberation  Geometry 

In  the  analyses,  there  was  a  requirement  to  develop  a  geometry  able  to  support 
both  monostatic  and  bistatic  reverberation  computations/measurements.  As  such,  a  16- 
element  VLA  was  used  in  the  modeling  effort.  Figure  4.1  shows  the  VLA  and  its 
geometry  relative  to  a  given  scattering  patch  at  horizontal  range  r.  We  see  that  the  mean 
horizontal  distance  of  the  VLA  to  a  small  scattering  patch  is  r,  in  see  Eq.  (2.43). 


Figure  4.1-  Geometry  of  VLA  and  Scattering  Patch 
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From  Figure  4.1,  we  see  that  the  16-element  VLA  spans  the  water  column  from  20  m  to 
80  m  with  the  array  elements  located  4  m  apart  in  depth.  A  single  element  at  48  m  depth 
(or  element  number  8  from  the  top)  was  chosen  as  the  source  with  all  16  elements 
receiving  the  reverberation. 

As  there  was  a  need  to  isolate  source/receiver  pairs,  MMPEREV  was  set  to 
generate  16  sets  of  output  files  using  each  array  element  as  a  reciprocal  source  and  then 
propagating  the  field  through  the  entire  space  with  the  environmental  and  source 
parameters  specified  in  the  next  sub-section. 


2.  The  Environment 

The  maximum  propagation  range  was  selected  to  be  5  km.  In  a  shallow  water 
environment,  the  mean  bottom  depth  was  assumed  to  be  100  m  with  no  mean  slope.  To 
perform  broadband  as  well  as  CW  analysis,  a  center  frequency  of  250  Hz  was  chosen 
with  a  250  Hz  bandwidth  divided  into  512  discrete  propagation  frequencies.  The 
frequency  span  was  from  124.7554  through  374.7554  Hz  computed  automatically  by 
MMPEREV.  Seven  control  and  environmental  input  files  were  required  by 
MMPE/MMPEREV  to  perform  the  computations.  In  order  to  properly  define  the 
environment  for  running  the  MMPEREV  model,  there  was  a  need  to  define  the  baseline 
features  of  frequency,  depth,  grid  (depth  and  range),  etc. 


In  order  that  the  time  steps  and  range  steps  coincide  appropriately  for  coherent 
addition  as  needed,  we  start  by  seeing  that  in  reduced  time  at  the  mth  range  step,  we  have 


Tm=tm  — —  , 


(4.1) 


where  Tm  is  reduced  time,  tm  is  true  time  at  range  step  m,  rm  is  the  range,  and  c0  is  the 
reference  sound  speed.  Hence, 


t  =7’+_2L- 


(4.2) 


Now,  true  time  at  rm+ 1  is 
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(4-3) 


^  ^,4.1  rr.  r  +  A  r 
tm+1  =T +-^  =  T +— - 


where  A r  is  the  range  step  size  (not  the  MMPEREV  computation  range  step  size,  which 
will  be  notated  as  Sr ),  which  leads  to 


^m+ 1 


The  frequency  step  is  defined  by 

A /=  — ,  (4.5) 

nf- 1 

where  nf  is  the  number  of  frequencies  we  would  like  to  compute  in  a  broadband  scenario 
with  BW  being  the  bandwidth  selected.  The  time  increment  then  becomes 

Al  =  — — —  =  j - —  =  (4.6) 

«/ -1  (n/-l)A/  BW 

Suppose  we  want  a  bandwidth  of  250  Hz,  a  center  frequency  of  250  Hz,  and  a 
frequency  step  size  of  about  V2  Hz,  then  selecting 


A  t  = - =  4  ms ,  (4.7) 

250 

we  get 

nf  =  512 ,  (4.8) 

250 

A/ =  —  =  0.49  Hz. 

511 

Note  that  nf  has  to  be  a  radix-2  integer  because  of  the  Fast-Fourier  Transform  (FFT)  that 
was  used  to  compute  spectral  components  in  MMPEREV  and  in  the  MATLAB  post¬ 
processing. 

We  now  constrain  the  range  step  to  be  an  integer  multiple  of  the  MMPEREV 
computational  range  step  size  by  defining 
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Ar  =  kxSr,  (4.9) 

where  it  is  an  arbitrary  integer  to  constrain  the  relationship  with  the  MMPEREV 
computation  range  step  size,  Sr .  By  forcing  this  constraint,  we  have  defined  the  spatial 
resolution  of  the  generated  solution  to  explicitly  overlap  the  temporal  resolution.  This 
will  allow  us  to  easily  combine  solutions  from  different  range  steps  coherently  in  the 
time-domain. 

To  resolve  time,  we  also  define 


—  =jxAt, 

c0 


(4.10) 


where  j  is  another  arbitrary  integer  to  constrain  the  time-range  relationship  for  the 
purpose  of  computations  in  MMPEREV.  For  At  =  4 ms,  selecting  j  =  k=  1  and 


c0  =1500 m/s,  we  get 


Ar  =  6m , 
Sr  =  6  m. 


(4.11) 


Hence  for  a  given  maximum  range  of  5  km,  the  number  of  range  steps  to  be  computed  is 


5000 

nr  = - 

6 


«  833 


(4.12) 


in  MMPEREV.  Note  that  the  maximum  computational  range  will  then  be 
=4.998  km. 


Finally,  we  select  the  number  of  total  computational  depth  points  to  be  256,  which 
would  also  determine  the  FFT-size  for  the  split-step  Fourier  algorithm,  also  a  radix-2 
integer.  This  provides  a  vertical  spatial  resolution  for  computation  of  roughly  3  m.  The 
following  table  summarizes  the  environmental  parameters  used  in  the  modeling  in 
MMPE,  noting  that  much  of  the  parameters  used  for  the  bottom  and  deep-bottom 
properties  were  assumed. 
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Filename/Parameter 

Value 

Remarks 

Main  Control  File:  pefiles.inp 

Number  of  depth  points 

256 

Radix-2  integer  required 

Minimum  depth 

0  m 

Maximum  depth 

Number  of  range  steps 

833 

Minimum  range 

0  m 

Maximum  range 

5.0  km 

Range  step  size 

6  m 

Maximum  computed  depth 

400  m 

Reference  sound  speed 

1500  m/s 

Source  depths 

Varying 

Array  elements  at  20,  24, 28,  32, 
36, 40, 44, 48, 52, 56, 60, 64, 

68, 72, 76  and  80  m  depths. 

Center  frequency 

250  Hz 

Frequency  bandwidth 

250  Hz 

No.  of  Frequencies 

512 

Radix-2  integer  required 

Sound  Speed  File:  pessp.inp 

Water  column  sound  speed 

1500  m/s 

Range  independent 

No.  ofSSPs 

1 

Bathymetry:  pebath.inp 

Mean  bottom  depth 

100  m 

Range  independent 

No.  of  depth  points 

1 

Bottom  properties:  pebotprop.inp 

Bottom  sound  speed 

1700  m/s 

Sound  speed  gradient 

0 

Relative  density 

1.0 

No  density  variations 

Compressional  attenuation 

0.2  dB/km/Hz 

Shear  speed 

0 

Not  modeled 

Shear  attenuation 

0 

Not  modeled 

Deep  Bottom  Bathymetry:  pedbath.inp 

Depth 

3000  m 

Deep  Bottom  Properties:  pedbotprop.inp 

Deep  bottom  sound  speed 

2000  m/s 

Sound  speed  gradient 

0 

Relative  density 

3.0 

No  density  variations 

Compressional  attenuation 

0.25  dB/km/Hz 

Shear  speed 

0 

Not  modeled 

Shear  attenuation 

0 

Not  modeled 

RMS  Perturbations  (input  to  MMPEREV 
during  run) 

Interface 

1  m 

As  explained  in  Chapter  III 

Volume 

15  m/s 

Table  4.1 -MMPEREV] 

iiput  Environmental  Properties 
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3.  MMPE  and  MMPEREV  Runs 

With  the  propagation  environment  specified,  MMPEREV  was  setup  to  make  16 
runs,  each  for  a  different  source  depth  from  20-80  m  at  4  m  increments  in  depth.  The 
ability  to  propagate  the  same  signal  out  from  the  various  receivers  to  the  scattering  patch 
is  based  on  reciprocity  as  explained  in  Chapter  II,  and  using  Eq.  (2.68b)  we  are  able  to 
arrive  at  the  reverberant  field.  As  such,  16  output  binary  files  for  the  configuration  above 
were  computed  and  stored.  In  order  to  make  comparisons  with  the  unperturbed  data, 
MMPE  was  also  setup  to  make  16  runs  for  the  same  array  element  depths  based  on  the 
same  environmental  conditions  shown  in  Table  4.1.  The  following  table  summarizes  the 
filenames  used. 


Receiver 

/Source 

Depth 

MMPEREV  Output  File 
(perturbed  data) 

MMPE  Output 
File  (unperturbed 
data) 

20m 

peoutzbrms  1  dcrms  1 5  s20.bin 

peouts20.bin 

24m 

peoutzbrms  1  dcrms  1 5  s24.bin 

peouts24.bin 

28m 

peoutzbrms  1  dcrms  15s28.bin 

peouts28.bin 

32m 

peoutzbrms  ldcrms  15s32.bin 

peouts32.bin 

36m 

peoutzbrms  Idcrms  15s36.bin 

peouts36.bin 

40m 

peoutzbrms  1  dcrms  1 5  s40.bin 

peouts40.bin 

44m 

peoutzbrms  Idcrms  15s44.bin 

peouts44.bin 

48m 

peoutzbrms  1  dcrms  1 5  s48  .bin 

peouts48.bin 

52m 

peoutzbrmsldcrmsl5s52.bin 

peouts52.bin 

56m 

peoutzbrms  1  dcrms  15s56.bin 

peouts56.bin 

60m 

peoutzbrms  1  dcrms  15s60.bin 

peouts60.bin 

64m 

peoutzbrms  1  dcrms  1 5s64.bin 

peouts64.bin 

68m 

peoutzbrms  1  dcrms  15s68  .bin 

peouts68.bin 

72m 

peoutzbrms  1  dcrms  15s72.bin 

peouts72.bin 

76m 

peoutzbrms  1  dcrms  1 5s76.bin 

peouts76.bin 

80m 

peoutzbrmsldcrmsl5s80.bin 

peouts80.bin 

Table  4.2  -MMPEREV/MMPE  Output  Files 


We  can,  from  a  single  CW  run  of  MMPE  and  MMPEREV,  make  qualitative 
comparisons  of  the  effects  of  the  interface  and  volume  perturbations. 
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Depth  (m)  Depth  (m) 


0.5  1  1.5  2  2.5  3  3.5  4  4.5  5 


Range  (km) 

Figure  4.2  -  Transmission  Loss  (dB  re  1  m)  with  No  Perturbation, 
Source:  48  m,  Frequency:  250  Hz 


0.5  1  1.5  2  2.5  3  3.5  4  4.5  5 


Range  (km) 

Figure  4.3  -  Transmission  Loss  (dB  re  1  m)  with  Perturbation, 
Source:  48  m,  Frequency:  250  Hz 
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When  the  perturbation  was  introduced,  the  structure  in  Figure  4.2  was  disrupted  as  shown 
in  Figure  4.3.  The  dispersed  penetration  into  the  volume  is  also  apparent.  In  the  next  few 
sections,  we  will  discuss  the  post-processing  implemented  to  compute  reverberation, 
vertical  correlation,  and  spectral  analysis. 

B.  POST-PROCESSING  1  -  REVERBERATION  (MATLAB 
IMPLEMENTATION) 

Recall  that  in  Chapter  II  we  derived  the  CW  reverberant  field  for  a  scattering 
patch  due  to  the  interface  and  volume  in  Eqs.  (2.33)  and  (2.61),  respectively.  In  order  to 
arrive  at  the  reverberant  field  at  the  receiver  due  to  a  co-located  or  displaced  source,  we 
had  to  find  the  interactions  of  these  fields  over  the  entire  scattering  area.  We  will  also 
perform  broadband  and  CW  computations  by  utilizing  the  same  output  from  MMPEREV 
set  to  run  on  the  various  geometrical  configurations  for  the  VLA  described  in  the 
previous  sub-section.  Again,  we  were  able  to  separate  the  bottom  interface  and  bottom 
volume  analyses  because  of  their  different  perturbations.  The  separate  analysis  also 
afforded  closer  discrimination  of  the  structure  of  one  independently  of  the  other. 

The  source  was  located  at  48  m  depth,  which  was  in  the  vicinity  of  the  center  of 
the  water  column.  Each  of  the  array  elements  (including  the  48  m  array  element)  were 
then  treated  as  independent  receivers.  Each  of  the  output  files  from  MMPEREV  and 
MMPE  could  be  paired  with  the  output  file  for  the  48  m  source  depth  binary  files  for 
computation.  Naturally,  the  48  m  source  and  48  m  receiver  files  were  the  same  binary 
file,  meaning  a  monostatic  geometry  while  the  other  transmitter/receiver  pairs  reflected  a 
vertical  bistatic  geometry. 

Both  perturbed  and  unperturbed  data  generated  using  MMPEREV  and  MMPE, 
respectively,  were  analyzed  to  provide  useful  comparisons  and  further  analyses.  We  will 
also  show  the  results  for  each  analysis  performed. 

1.  Two-way  Mono/Bistatic  CW  Reverberation  Analysis 

For  the  interface,  we  have  from  Eq.  (2.33)  and  subsequently  in  Eq.  (2.43), 
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(4.13) 


(|P-f  )  =p-K(r)|Va(r)|Jt>%,20»)AA.  ■ 

This  is  simply  proportional  to  the  magnitude-squared  of  the  forward  and  backward 
propagated  field  functions.  As  the  perturbation  spectrum  and  other  field  constants  are 
generated  by  MMPEREV,  the  MATLAB  post-processing  that  computes  the  CW 
reverberation  is  just  the  multiplication  of  the  magnitude-squared  field  functions  with 
cylindrical  spreading. 


Wn  (Ol^Oa  (r) 

Vr  1  4r 

and  the  interface  reverberation  loss  (RLb)  is  then 


(4.14) 


RLb(r)  =  A -20 log 

r 

where  A  is  a  constant  which  accounts  for  the  other  parameters  in  Eq.  (2.50). 


(4.15) 


The  volume  was  also  treated  as  described  in  Eq.  (2.66b)  and  the  volume 
reverberation  loss  may  be  written  as 


RLv(r)  =  B-20log 


-  \n{r,z)\yfT(.r,z)\yR(,r,z)\dz 

^  z>zb 


(4.16) 


where  B  is  another  constant  which  accounts  for  the  other  terms  in  Eq.  (2.66b).  Using  now 
the  same  16  sets  of  results  from  a  single  MMPEREV  and  another  16  sets  of  results  from 
the  MMPE  runs,  we  post-processed  the  field  function  information  to  arrive  at  the  CW 
interface  and  volume  analyses.  We  have  selected  the  minimum,  center  and  maximum 
frequencies  for  analysis.  The  MATLAB  files  written  to  post-process  the 
MMPEREV/MMPE  output  files  to  arrive  at  the  CW  analysis  are  given  in  Table  4.3.  The 
output  from  each  run  of  peout2way_bistat_r.m  was  saved  as  a  MATLAB  data  file  so  that 
the  data  set  could  be  loaded  into  the  MATLAB  workspace  without  having  to  re-process 
the  MMPEREV/MMPE  output  binary  files.  The  filenames  saved  are  listed  in  Table  4.5. 
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MATLAB  Post-Processing  Files 
Description 

Perturbed  Data 
Processing  Files 

Unperturbed  Data 
Processing  Files 

Binary  file  data  extractor:  Opens  both  source  and 
receiver  binary  files  (output  from  MMPE/ 
MMPEREV)  and  all  modeling  data.  Stores 
them  on  the  MATLAB  workspace. 

peoutl_bistat_r.m 

peout  l_bistat_np_r.m 

Interface  and  Volume  Reverberation  loss 
Computation:  Extracts  field  data  source  and 
receiver  files  at  selected  frequencies. 

Computes  the  two-way  reverberation  loss  for 
the  interface  and  volume  as  described  above. 
Extracts  field  data  at/mfn,  /r^r,  and  /™*. 

peout2way_bistat_r.m 

peout2  way_bi  stat„np_r .  m 

Table  4.4  -  MATLAB  Filenames  Created  for  CW  Reverberation  Analysis 


Receiver 

Depth 

Interface 
Perturbed  Data 
Ouput 

Volume 
Perturbed  Data 
Ouput 

Interface 

Unperturbed  Data 
Ouput 

Volume 

Unperturbed  Data 
Ouput 

20m 

bistat4820int_r.mat 

bistat4820int_np_r.mat 

bistat4820vol_np_r.mat 

24m 

bistat4824int_r.mat 

bistat4824voLr.mat 

bistat4824int_np_r.mat 

bistat4824vol_np_r.mat 

28m 

bistat4828int_r.mat 

bistat4828vol_r.mat 

bistat4828int_np_r.mat 

bistat4828vol_np_r.mat 

32m 

bistat4832int_r.mat 

bistat4832int_np_r.mat 

bistat4832vol_np_r.mat 

36m 

bistat4836int_r.mat 

bistat4836int__np_r.mat 

bi  stat48 3 6 vol_np_r.  mat 

40m 

bistat4840int_r.mat 

bistat4840vol_r.mat 

bistat4840int_np_r.mat 

bistat4840voLnp_r.mat 

44m 

bistat4844int_r.mat 

■JUMMUUM— 

bistat4844int_np_r.mat 

bistat4844vol_np_r.mat 

48m 

bistat4848int_r.mat 

bistat4848int_np_r.mat 

bistat4848vol_np_r.mat 

52m 

bistat4852int_j\mat 

bistat4852vol_r.mat 

bistat4852int_np_r.mat 

bistat4852vol_np_r.mat 

56m 

bistat4856int_r.mat 

bistat4856int_np_r.mat 

bistat4856vol_np_r.mat 

60m 

unamiiuH 

bistat4860int_np_r.mat 

bistat4860vol_np_r.mat 

64m 

Kgmwssm 

bistat4864int_np„r.mat 

bistat4864vol_np_r.mat 

bistat4868int_r.mat 

bistat4868int_np_r.mat 

bistat4868vol_np_r.mat 

■E35HI 

bistat4872int_r.mat 

bistat4 8 72 int_np_r .  mat 

bistat4872voI_np_r.mat 

76m 

bistat4876int_r.mat 

^  i  j  ^  Mi  m  k  I'M 

bistat4876int_np_r.mat 

bistat487 6vol_np_r.mat 

80m 

bistat48 80int_r.mat 

bistat4880int_np_r.mat 

bistat4880vol_np_r.mat 

Table  4.5  -  Processed  Output  Data  Files  in  CW  Analysis 


2.  CW  Analysis  Results 

A  typical  CW  reverberation  loss  plot  shows  a  20  log  r  (two-way  cylindrical 
spreading)  drop  in  range.  Figures  4.4  and  4.5  show  the  interface  and  volume 
reverberation  loss  for  a  source  depth  of  48  m,  receiver  depth  of  40  m  and  frequency  of 
250  Hz.  In  all  cases,  the  constant  scaling  parameters  A  and  B  were  set  equal  to  zero.  With 
the  same  source  and  receiver  depths  and  at  the  same  center  frequency,  we  were  able  to 
observe  finer  structures  in  the  interface  reverberation  (Figure  4.4)  compared  to  the 
volume  reverberation  plot  (Figure  4.5).  This  was  because  multiple  incoherent 
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multiplications  were  performed  at  all  depths  and  then  summed  vertically  m  order  to 
arrive  at  the  mean  squared  RPL  at  the  receiver  due  to  the  volume.  Note  that  all  data  sets 
show  very  little  energy  interacting  with  the  interface  or  volume  at  very  short  ranges,  i.e. 
<  0.1km .  This  is  due  to  the  wide-angle  limitation  of  the  parabolic  approximation. 

Source  depth:48m,  Receiver  depth:40m,  Frequency:  250Hz 

70  j— - 1  I  I  I  I  I  .  i .  HZ  i  . . 

LJi  —  Interface  reverberation  loss 

in  —  20log  r  drop  off 


Range  (km) 


Figure  4.4-  Interface  Reverberation  Loss  for  Two-Way  Transmission  With  Perturbation 


Source  depth:48m,  Receiver  depth:40m,  Frequency:  250Hz 


Interface  reverberation  loss 
20log  r  drop  off  _ 


Figure  4.5  -  Volume  Reverberation  Loss  for  Two-Way  Transmission  With  Perturbation 


A  comparison  between  perturbed  and  unperturbed  data  was  also  done.  Figures  4.6 
and  4.7  show  the  difference  (perturbed  RL  minus  unperturbed  RL)  to  see  the  fine-scale 
structures  due  to  the  introduction  of  the  perturbation.  Apparent  in  these  plots  are 
repeating  structures  at  about  every  0.5  km  for  ranges  beyond  about  2  km. 

(Perturbed  -  Unperturbed)  Source  Depth:  48m,  Receiver  depth:  40m,  Frequency:  250Hz 


Figure  4.6  -  Interface  Reverberation  Loss  Difference  Vs  Range 
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(Perturbed  -  Unperturbed)  Source  depth:  48m,  Receiver  depth:  40m,  Frequency:  250Hz 


Figure  4.7  -  Volume  Reverberation  Loss  Difference  Vs  Range 


Figures  4.8  and  4.9  show  the  interface  and  volume  reverberation  loss, 
respectively,  across  depths  plotted  against  range  for  the  center  frequency  of  250  Hz.  It 
was  apparent  that  the  interface  reverberation  showed  more  structure  while  volume 
reverberation  seemed  more  ‘blotchy’  and  revealed  little  structure.  This  smudging  of  data 
was  manifest  of  the  multiple  incoherent  (magnitude)  multiplications  computed  before 
vertically  summing  the  reverberation  in  depth.  What  was  apparent  however,  was  the 
‘201og  r'  drop  in  range  in  both  plots. 
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Range  (km) 


Figure  4.8  -  Interface  Reverberation  Loss  (dB  re  1  m) 
Across  Depth  Vs  Range 

Source  Depth:  48m,  Frequency:  250Hz 


1.5  2  2.5  3 

Range  (km) 


Figure  4.9  -  Volume  Reverberation  Loss  (dB  re  1  m) 
Across  Depth  Vs  Range 


3.  Two-way  Mono/Bistatic  Time-Domain  Reverberation  Analysis 

Time-domain  analyses  for  both  the  interface  and  volume  reverberation  were 
explained  in  Chapter  II,  Section  C.  In  essence,  multiple  frequencies  were  used  to  simulate 
a  pulsed  propagation.  We  will  use  Eqs.  (2.68b)  and  (2.68c)  to  formulate  the  MATLAB 
implementation  for  the  interface  and  Eqs.  (2.70a)  to  (2.72)  for  the  volume.  MMPE  was 
also  run  based  on  the  same  16  setups  for  the  VLA  to  provide  comparison  data  for 
unperturbed  data.  Implementation  was  completed  using  MATLAB  with  the  files  created 


as  follows: 


MATLAB  Post-Processing  Files 
Description 

Perturbed  Data 
Processing  Files 

Unperturbed  Data 
Processing  Files 

Binary  file  data  extractor:  Opens  both  source 
and  receiver  binary  files  (output  from 
MMPE/MMPEREV).  Extracts  all  modeling 
parameters  and  stores  them  on  the 

MATLAB  workspace. 

Peoutl_bistat.m 

peoutl_bistat_np.m 

Interface  and  Volume  Reverberation 

Computation:  Extracts  from  binary  files  the 
field  function  data  for  source  and  receiver 
and  computes  the  two-way  reverberation 
loss  for  die  interface  and  volume  as 
described  above. 

Peout2way_J)istat.m 

peout2way_bistat_np.m 

Table  4.6  -  MATLAB  Filenames  Created  for  Time-Domain  Reverberation  Analysis 

The  output  from  each  run  of  peout2way_bistat.m  was  saved  as  a  MATLAB  data  file  so 
that  the  data  set  could  be  loaded  into  the  MATLAB  workspace  without  having  to  re¬ 
process  the  MMPEREV/MMPE  output  binary  files.  These  are  listed  in  Table  4.7. 

Additionally,  peoutl_bistat.m  and  peout2way_bistat.m  were  modified  into 
function  form  in  order  to  allow  unattended  batch  processing  of  the  transmitter/receiver 
pairs.  MATLAB  routines  created  were:  master_run.m,  peoutl_bifn.m, 
peout2way_biintfn.m,  and  peout2way_bivolfn.m,  where  the  interface  and  volume 
processing  in  the  original  peout2way_bistat.m  was  split  up  into  two  files  for  processing 
efficiency. 
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Receiver 

Depth 

Interface 
Perturbed  Data 
Ouput 

Volume 
Perturbed  Data 
Ouput 

Interface 
Unperturbed 
Data  Ouput 

Volume 
Unperturbed 
Data  Ouput 

20m 

bistat4820int.mat 

bistat4820vol.mat 

bistat4820int_np.mat 

bistat4820voLnp.mat 

24m 

bistat4824int.mat 

bistat4824vol.mat 

bistat4824int_np.mat 

bistat4824vol_np.mat 

28m 

bistat4828int.mat 

bistat4828vol.mat 

bistat4828int_np  .mat 

bistat4828vol_np.mat 

32m 

bistat4832intmat 

bistat4832vol.mat 

bistat4832int_np.mat 

bistat4832vol_np.mat 

36m 

bistat4836int.mat 

bistat4836vol.mat 

bistat4836int_np.mat 

bistat4836vol_np.mat 

40m 

bistat4840int.mat 

bistat4840vol.mat 

bistat4840int_np.mat 

bistat4840vol_np.mat 

44m 

bistat4844int.mat 

bistat4844vol.mat 

bistat4844int_np.mat 

bistat4844vol_np.mat 

48m 

bistat4848int.mat 

bistat4848vol.mat 

bistat4848int_np.mat 

bistat4848vol_np.mat 

52m 

bistat4852int.mat 

bistat4852vol.mat 

bistat4852int_np.mat 

bistat4852vol_np.mat 

56m 

bistat4856int.mat 

bistat4856vol.mat 

bistat4856int_np.mat 

bistat4856vol_np.mat 

60m 

bistat4860int.mat 

bistat4860vol.mat 

bistat4860int_np.mat 

bistat4860vol_np.mat 

64m 

bistat4864int.mat 

bistat4864vol.mat 

bistat4864int_np.mat 

bistat4864vol_np.mat 

68m 

bistat4868int.mat 

bistat4868vol.mat 

bistat4868int_np.mat 

bistat4868vol_np.mat 

wbbfr m 

bistat4872int.mat 

bistat4872vol.mat 

bistat4872int_np.mat 

bistat4872vol_np.mat 

bistat4876int.mat 

bistat4876vol.mat 

bistat4876int_np.mat 

bistat4876vol_np.mat 

bistat4880int.mat 

bistat4880vol.mat 

bistat4880int_np.mat 

bistat4880vol_np.mat 

Table  4.7  -  Processed  Output  Data  Files  in  Time-Domain  Analysis 


4.  Time-Domain  Analysis  Results 

Having  discussed  the  implementations,  we  now  show  some  results  obtained.  First, 
we  looked  at  the  general  structure  of  interface  and  volume  reverberation  in  the  time- 
domain  at  a  particular  depth  chosen  arbitrarily  as  40  m.  Showing  very  much  the  same 
characteristics,  Figures  4.10  and  4.11  depict  the  time-domain  interface  and  volume 
reverberation  loss,  respectively.  The  figures  provide  an  understanding  of  the  structure  of 
the  reverberation  return  of  a  pulsed  transmission  with  returns  lasting  for  the  entire  extent 
of  the  calculations  both  in  the  interface  and  volume. 
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Figure  4.10  -  Time-Domain  Interface  Reverberation  Loss 


Time-domain  data  sets  were  also  used  to  provide  useful  comparisons  with 
unperturbed  data  sets.  To  do  so,  the  time  axis  was  ‘converted’  to  a  range  axis  by  use  of 
the  reference  sound  speed,  c0 ,  such  that  range  =  c0  x  time .  This  produced  a  range  axis 

from  0  to  5  km.  The  data  comparison  then  proceeded  along  the  same  path  as  that  in  the 
CW  analysis  in  the  previous  section.  We  first  show  sample  time-series  computed 
reverberation-range  plots  for  the  interface  and  volume.  Shown  in  Figures  4.12  and  4.13 
are  the  time-domain  reverberation  loss  plots  equivalent  to  Figures  4.4  and  4.5  for  the  CW 
analyses.  Again,  the  ‘201og  r’  drop  off  was  apparent. 


Source  depth:  48m,  Receiver  depth:  40m 


Figure  4.12  -Time-Domain  Converted  Interface  Reverberation  Loss 
(dB  re  1  m)  Vs  Range,  With  Perturbation 
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Source  depth:  48m,  Receiver  depth:  48m 


Figure  4.13  -  Time-Domain  (Converted  to  Range) 
Volume  Reverberation  Loss  Vs  Range,  With  Perturbation 


Now  looking  across  depths,  we  observe  multipath  structures  shown  in  Figures 
4.14  and  4.15  due  to  both  interface  and  volume  perturbations,  as  expected.  Spatial 
disposition  of  each  source  was  important  in  observing  the  arrival  time  of  the  first  returns, 
most  apparent  at  the  receiver  element  at  20  m  depth,  while  arrivals  were  earlier  at  the 
receiver  element  at  80  m  depth.  This  was  due  simply  to  the  closer  proximity  of  the  deep 
receiver  to  the  bottom  thereby  receiving  returns  earlier  than  the  shallow  receivers.  No 
other  structures  were  apparent  from  these  two  figures. 
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Source  Depth:  48m 


Time  (sec) 


Figure  4.14  -  Interface  Reverberation  Loss  Across  Depth  Vs  Time 


Source  Depth:  48m 


Time  (sec) 


Figure  4.15  -  Volume  Reverberation  Loss  Across  Depth  Vs  Time 
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5.  MATLAB  Data  Analysis 

As  a  final  note  on  the  work  done  in  both  data  and  signal  analyses,  several 
MATLAB  files  were  also  written  to  ensure  that  the  data  analysis  was  complete.  These  are 
found  in  Table  4.8.  In  order  to  derive  further  value,  vertical  correlations  and  spectral 
analyses  were  performed.  These  are  discussed  in  the  next  sections. 


MATLAB  Post-Processing  Files  Description 

MATLAB  Files 

Range-Series  Data  Processor:  Using  data  obtained  from 
the  output  file  created  by  peout31_pre_bistat_r.m  and 
peout31_pre_bistat_r_np.m  respectively,  this  routine 
computes  and  displays  reverberation  loss  and 
compares  these  RLs  between  perturbed  and 
unperturbed  range-series  data. 

range_series_datanalysis.  m 

Time-Series  Data  Processor:  Using  data  obtained  from 
the  output  file  created  by  peout31_pre_bistat_r.m  and 
peout31_pre_bistat_r_np.m  respectively,  this  routine 
computes  and  displays  reverberation  loss  and 
compares  these  RLs  between  perturbed  and 
unperturbed  range-series  data. 

time_series_datanalysis.m 

Table  4.8  -  MATLAB  Routines  Created  for  Data  Analyses  and  Comparisons 


C.  POST-PROCESSING  2  -  VERTICAL  CORRELATION  &  PEAK 
CORRELATION  IN  CW  AND  TIME-DOMAIN 

Having  completed  CW  and  time-domain  analyses,  correlation  between  the 
monostatic  (source  and  receiver  co-located  at  48  m)  reverberation  data  set  with 
reverberation  data  sets  from  all  other  receiver  depths  was  performed.  The  data  files  saved 
from  the  CW  (Table  4.5)  and  time-domain  (Table  4.7)  analyses  were  used  to  compute  the 
vertical  correlations  for  both  CW  (range  correlation)  and  time-domain  (temporal 
correlation).  MATLAB  routines  were  written  to  perform  the  vertical  correlations.  All 
vertical  correlations  are  computed  relative  to  a  source  depth  of  48  m.  Thus,  the 
autocorrelation  would  occur  at  the  receiver  depth  of  48  m,  giving  the  theoretical 
maximum  correlation  over  depth.  These  were  done  for  both  the  interface  and  volume 
with  and  without  perturbations.  Selected  results  are  shown  in  this  section. 
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1.  Vertical  Correlation  in  Range  &  Peak  Vertical  Correlation 

We  first  looked  at  CW  vertical  correlation.  Table  4.9  lists  the  MATLAB  files 
created  to  compute  the  vertical  correlations.  We  have  selected  a  source  at  48  m  and  at  the 
center  frequency  of  250  Hz  for  this  analysis.  With  the  perturbation  included,  Figure  4.16 
shows  the  vertical  correlation  for  the  interface  while  Figure  4.17  shows  that  in  the  case  of 
the  volume.  It  should  be  noted  that  due  to  the  incoherent  processing  of  the  two-way 
pressure  fields,  the  first  0.5  km  does  hot  accurately  reflect  the  true  reverberation  loss. 
Hence,  the  first  0.5  km  has  been  truncated  prior  to  performing  the  volume  vertical 


correlations  shown  in  Figure  4.17. 


MATLAB  Post-Processing  Files 
Description 

MATLAB  Files 
(Perturbed) 

MATLAB  Files 
(Unperturbed) 

Correlation  Computation:  Extracted  data  from 
Table  4.5  and  computed  vertical  and  peak 
vertical  correlations  with  the  data  at  source 
depth,  48m. 

Peout3_bistat_r.m 

peout3_bistat_r_np  .m 

Data  Processing:  Extracted  data  from  Table  4.5 
and  stored  workspace  variables  into  an 
output  data  file  for  future  analyses. 

Peout3  l_pre_bistat_r.m 

peout3 1  _pre_bistat_r_np.m 

Correlation  Computation:  Computes  vertical 
and  peak  vertical  correlations  as  in 
peout3_bistat_r.m  and 
peout3_bistat_r_np.m  respectively,  but  uses 
output  file  created  by 
peout31_pre_bistat_r.m  and 
peout3  l_pre_bistat_r_np.m  respectively. 

Peout3  l_bistat  jr.m 

peout3  l_bistat_r_np.m 

Output  filename  created  by 
peout3  l_pre__bistat_r.m  and 
peout3  l_pre_bistat_r_np.m  respectively 

Bistat_depmg_data.mat 

bistat_depmg_np_data.mat 

Table  4.9  -  MATLAB  Filenames  Created  for  Computing  Vertical  Correlations 
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Relative  Depth  (m)  Relative  Depth  (m) 


Source  Depth:  48m,  Frequency:  250Hz 


Range  Lag  (km) 

Figure  4. 16  -  Vertical  Correlation  of  Interface  Reverberation  Loss 
in  Range  Vs  Relative  Depth 


Source  Depth:  48m,  Frequency:  250Hz 


-3-2-10123 

Range  Lag  (km) 


Figure  4.17  -  Vertical  Correlation  of  Volume  Reverberation  Loss 
in  Range  Vs  Relative  Depth  (first  0.5 km  truncated) 
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It  was  apparent  that  finer  scale  structures  were  present  in  the  interface,  but  the  volume 
showed  strong  correlation  throughout  depth  for  much  of  the  computations  in  range  up  to 
about  3  km  due  to  the  smoothing  effect  of  the  incoherent  processing  of  the  CW  signal. 

Peak  correlation  values  were  also  extracted  from  the  above  vertical  correlation 
structures  in  order  to  see  how  rapidly  the  signal  decorrelated  over  depth.  As  would  be 
expected  from  the  above  analysis,  the  volume  would  tend  to  stay  more  correlated  over 
depth  while  the  interface  showed  more  decorrelation  over  depth.  However,  even  the 
interface  reverberation  loss  remained  quite  correlated  with  minimum  correlation  values 
down  only  about  1  dB.  This  is  also  due  presumably  to  the  incoherent  processing  involved 
in  the  derivation  of  the  mean  squared  RPL.  Figures  4.18  and  4.19  show  a  comparison  of 
the  peak  vertical  correlations  between  perturbed  and  unperturbed  data.  While  it  is 
apparent  that  the  volume  reverberation  of  the  CW  signal  tends  to  decorrelate  slower  than 
the  interface  in  both  scenarios,  the  addition  of  the  perturbations  does  not  appear  to  affect 
the  vertical  correlations  significantly. 


Source  Depth:  48m,  Frequency:  250Hz 

-30 |  i  i  i  i  i  i  i  r 


-1.6  -1.4  -1.2  -1  -0.8  -0.6  -0.4  -0.2  0  0.2 

Peak  Correlation  (dB) 

Figure  4.18  -  Peak  Vertical  Correlations  (dB)  of 
Perturbed  Data  Over  Relative  Depth 
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-30 


Source  Depth:  48m,  Frequency:  250Hz 


i - 1 - : - 1 - 1 - 1 - 1 - r 


40 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 

-1.6  -1.4  -1.2  -1  -0.8  -0.6  -0.4  -0.2  0  0.2 

Peak  Correlation  (dB) 

Figure  4.19  -  Peak  Vertical  Correlations  (dB)  of 
Unperturbed  Data  Over  Relative  Depth 

2.  Vertical  Temporal  Correlation  &  Peak  Vertical  Temporal 
Correlation 

As  with  the  CW  signals,  MATLAB  files  were  created  to  compute  the  vertical 
correlations  in  the  time-domain  for  the  broadband  signals.  These  are  listed  in  Table  4.9. 
From  the  results  of  the  vertical  temporal  correlations  and  peak  vertical  temporal 
correlations,  we  obtain  Figures  4.20  through  4.23. 

Apparent  from  Figures  4.20  and  4.21  are  finer  scale  structures  seen  in  both  the 
interface  and  volume  reverberation  level.  Of  importance  was  the  fact  that  the  coherent 
convolution  done  in  arriving  at  the  time-domain  analysis  preserved  the  coherent  time 
structures,  in  contrast  to  the  CW  analyses. 
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MATLAB  Post-Processing  Files 
Description 

MATLAB  Files 
(Perturbed) 

MATLAB  Files 
(Unperturbed) 

Correlation  Computation:  Extracted  data  from 
Table  4.7  and  computed  vertical  and  peak 
vertical  correlations  with  the  data  at  source 
depth,  48m. 

peout3_bistat.m 

peout3_bistat_np.m 

Data  Processing:  Extracted  data  from  Table  4.7 
and  stored  workspace  variables  into  an 
output  data  file  for  future  analyses. 

peout3  l_pre_bistat.m 

peout3  l_pre_bistat_np.m 

Correlation  Computation:  Computes  vertical 
and  peak  vertical  correlations  as  in 
peout3_bistat.m  and  peout3_bistat_np.m 
respectively,  but  uses  output  file  created  by 
peout3  l_pre_bistat.m  and 
peout3  l_prej)istat_np.m  respectively. 

peout31_bistat.m 

peout3  l_bistat_np.m 

Output  filename  created  by 
peout31_pre_bistat.m  and 
peout3  l_pre_bistat_np.m  respectively 

bistat_deptime_data.mat 

i 

bistat_deptime_np  _data.mat 

Table  4.10  -  MATLAB  Filenames  Created  for  Computing  Vertical  Correlations 


Figure  4.20  -  Vertical  Temporal  Correlation  of 
Interface  Reverberation  Loss  Vs  Relative  Depth 
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Source  Depth:  48m 


Time  Lag  (sec) 

Figure  4.21  -  Vertical  Temporal  Correlation  of 
Volume  Reverberation  Loss  Vs  Relative  Depth 

Peak  vertical  temporal  correlations  were  extracted  for  the  interface  and  volume 
relative  to  a  48m  source,  and  the  results  are  displayed  in  Figure  4.22.  Similar  analysis 
was  performed  on  data  computed  without  the  perturbations  and  these  results  are 
presented  in  Figure  4.23.  With  the  same  horizontal  scales,  Figures  4.22  and  4.23  could  be 
compared  to  see  the  effects  of  introducing  perturbations.  The  perturbations  appear  to 
affect  the  volume  reverberation  loss  vertical  correlation  more  than  the  interface 
reverberation  loss,  although  the  effect  is  minor  in  both  cases.  For  both  perturbed  and 
unperturbed  environments,  the  vertical  3dB  decorrelation  length  scale  is  on  the  order  of  a 
few  meters. 
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Figure  4.22  -  Peak  Vertical  Temporal  Correlations 
of  Perturbed  Data 


Source  Depth:  48m 


Figure  4.23  -  Peak  Vertical  Temporal  Correlations 
of  Unperturbed  Data 
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D.  POST-PROCESSING  3  -  SPECTRAL  ANALYSES 


In  this  final  section,  we  investigate  the  spectral  content  of  the  reverberation 
pressure  obtained  for  both  the  CW  and  broadband  signals.  There  was  a  need  to  convert 
the  time-domain  into  range  using  the  reference  sound  speed,  Co,  like  the  conversion  done 
in  Section  3.B  in  the  time-domain  analysis.  As  before,  comparison  with  unperturbed  data 
was  also  done.  MATLAB  routines  were  again  written  and  tabulated  as  follows. 


MATLAB  Post-Processing  Files 
Description 

MATLAB  Files 

Computes  the  Fourier  Transform  of  the  CW 
data  obtained  from  the  output  file  created  by 
peout3  l_pre_bistat.m  and 
peout31_pre_bistat_np.m  respectively.  This 
analysis  has  built  within  itself  a  comparison 
with  unperturbed  data  already. 

range_series_signalysis.m 

Computes  the  Fourier  Transform  of  the 

Converted  Time-Domain  data  obtained  from 
the  output  file  created  by 
peout31_pre_bistat.m  and 
peout31_pre__bistat_np.m  respectively.  This 
analysis  has  built  within  itself,  a  comparison 
with  unperturbed  data  already. 

time_series_signalysis.m 

Table  4.11  -  MATLAB  Routines  Created  for  Spectral  Decomposition 


Assuming  a  monostatic  geometry,  we  recall  from  the  analyses  in  Chapter  II  and 
from  Eqs.  (4.15)  and  (4.16)  that  the  interface  reverberation  loss  may  be  defined  as 

RI^  =  2 TLb  +  constants  ,  (4- 17) 

where 


TLb  =  10  log|p_  (r)|2  =  20  log|^  (r)|  - 10  log  r 
since,  by  definition, 

\p(r)\  =  -j=\y(r)\. 

Hence 

RLh  °c  20  log|^  (r)|2  -  20  log  r . 

We  define  the  range-reduced  reverberation  loss  as 


(4.18) 


(4.19) 


(4.20) 
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(4.21) 


RLbjeiuced  =RLb  +  20 log  r  oc  20  log|^  (r)|2 . 

Now  we  assume  that  we  can  represent  the  field  function  in  terms  of  a  spectrum, 
i.e. 

\vM'*\W(K)e‘K'dK.  (4'22) 

It  is  possible  then  that  the  spectral  content  of  the  field  can  be  related  to  the  statistical 
characteristics  of  the  perturbations  and  can  thus  be  extracted  with  the  use  of  some  signal 
processing  tools.  In  other  words,  the  spectrum  W(K)  may  be  related  to  the  perturbation 
spectra. 

The  main  approach  was  the  use  of  discrete  Fourier  transforms  (DFT)  in  an 
attempt  to  extract  the  spectral  components  of  the  reverberation  data.  They  are: 

•  Signal  Analysis  1:  Power  Spectral  Density  -  DFT  of  the  magnitude  squared 
range-reduced  reverberation  (i.e.  r2x|p_|2). 


•  Signal  Analysis  2:  Power  Ratio  Spectral  Density  -  DFT  of  the  magnitude 
squared  of  the  ratio  of  the  reverberation  (perturbed  divided  by  unperturbed), 


i.e., 


\P~  perturbed 


P -unperturbed 


2  * 


Recall  that  the  typical  reverberation  loss  observed  from  the  range  series  (CW) 
data  for  the  interface  and  volume  were  presented  earlier  in  this  chapter  in  Figures  4.4  and 

4.5  respectively.  The  time-domain  reverberation  loss  data  was  correspondingly  shown  in 
Figures  4.10  and  4.11.  The  difference  in  RL  (ratio)  between  the  perturbed  and 
unperturbed  data  (CW  analysis)  for  the  interface  and  volume  were  also  shown  in  Figures 

4.6  and  4.7.  Figures  4.24  and  4.25  show  the  CW  reverberation  loss  while  Figures  4.26 
and  4.27  show  the  time-domain  (converted  to  range)  reverberation  loss. 
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Source  depth:  48,  Receiver  depth:  40m,  Frequency:  250Hz 


Figure  4.24  -  Interface  Reverberation  Loss  vs 
Range  Using  Range  Series  Data 


Source  depth:  48,  Receiver  depth:  40m,  Frequency:  250Hz 


Figure  4.25  -  Volume  Reverberation  Loss  vs 
Range  Using  Range  Series  Data 

When  using  time-series  data  to  perform  this  analysis,  time  was  converted  to  range 
in  order  to  relate  length  scales  of  the  reverberation  with  perturbation  scales.  It  is  perhaps 
interesting  to  note  that  the  time-series  data  showed  no  appreciable  differences  between 
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the  perturbed  and  unperturbed  data  sets,  in  contrast  to  the  CW  range-series  data  displayed 
above.  The  reason  for  this  is  currently  unknown. 


(Time  Series)  Source  depth:  48m,  Receiver  depth:  40m 


Figure  4.26  -  Interface  Reverberation  Loss  (dB  re  1  m)  vs 
Range  Using  Time  Series  Data 

(Time  Series)  Source  depth:  48m,  Receiver  depth:  40m 


— 

— 

- 1 - 

- r 

—  Unperturbed  reverberation  loss 

—  Perturbed  reverberation  loss 

1.  Signal  Analysis  1  -  Power  Spectral  Density 

The  magnitude-squared  of  the  range-reduced  perturbed  reverberation  was 
analyzed  using  the  DFT,  giving 

PSD  =  DFT{\p_\2  x  r2 }  (4.23) 

where  p-  is  the  reverberation  field  of  the  interface  or  volume.  Figures  4.28  and  4.29  show 
the  normalized  power  spectrum  of  the  range-reduced  interface  and  volume 
reverberations,  respectively.  Note  that  these  plots  are  presented  with  the  wavenumber 
axis  displayed  on  a  logarithm  scale  in  order  to  elucidate  exponential  behavior.  While 
structurally  the  same,  i.e.,  lower  wavenumbers  predominating  the  spectra,  strong  fine 
scale  structures  were  evident  in  both  plots  for  the  perturbed  data.  This  alludes  to  the 
presence  of  perturbations  which  may  be  extracted  from  the  spectra.  In  the  volume  plot, 
Figure  4.29,  the  perturbed  data  distinctly  showed  less  drop-off  beyond  a  wavenumber  of 
roughly  0.2  m'1  corresponding  to  wavelengths  smaller  than  approximately  30  m. 


Source  depth:  48m,  Receiver  depth:  40m,  Frequency:  250Hz 


Wavenumber,  (nrf1) 


Figure  4.28  -  Normalized  Power  Spectrum  of 
Range-Reduced  Interface  Reverberation  (Range-Series) 
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Source  depth:  48m,  Receiver  depth:  40m,  Frequency:  250Hz 


Wavenumber,  (nrf1) 

Figure  4.29  -  Normalized  Power  Spectrum  of 
Range-Reduced  Volume  Reverberation  (Range-Series) 

Figures  4.30  and  4.31  show  the  results  for  the  broadband  signal  (after  converting 
time  to  range)  for  the  interface  and  volume.  Both  fine  and  large  scale  structures  are 
observed.  While  repeating  structures  were  mostly  preserved  when  the  perturbation  was 
included,  the  longer  scale  perturbations  may  have  augmented  the  spectra. 
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2.  Signal  Analysis  2  -  Power  Ratio  Spectral  Density 

The  quotient  of  the  perturbed  over  the  unperturbed  reverberation  data  was 
magnitude-squared  to  provide  a  comparison  of  the  levels.  This  was  then  Fourier 
transformed  in  order  to  reveal  any  special  spectral  content.  The  power  ratio  spectral 
density  (. PRSD )  was  then  defined  as 

2  ’ 

PRSD  =  DFT |  -perturbed  2  1.  (4-24) 

P -unperturbed 

Figures  4.32  and  4.33  show  the  spectral  components  of  the  ratio  between  the  perturbed 
and  unperturbed  data  sets  for  the  CW  signal.  While  it  was  difficult  to  make  any  apparent 
observations,  the  interface  reverberation  data  showed  a  rather  flat  spectrum  while  the 
volume  reverberation  data  showed  a  roll-off  at  higher  wavenumbers.  The  physical 
justification  for  this  is,  as  yet,  unclear. 


(IPerturbedf/IUnperturbedl2)  Source  depth:  48m,  Receiver  depth:  40m,  Frequency:  250Hz 


(|Perturbed|2/|Unperturbed|2)  Source  depth:  48m,  Receiver  depth:  40m,  Frequency:  250Hz 
Oi - 1 - 1 — * — ;  rrm - : - t — : — i  rm - : - : — ! — " 


Wavenumber,  (m_1) 


Figure  4.33  -  Normalized  Spectrum  of  the  Volume  Power  Ratio 
Using  Range-Series  Data 

Finally,  we  computed  the  power  ratio  spectral  density  based  on  the  time-series 
data,  as  shown  in  Figures  4.34  and  4.35  for  the  interface  and  volume,  respectively.  These 
plots  were  consistent  with  the  above  CW  analysis  for  the  ratio  of  the  perturbed  to  the 
unperturbed  data.  The  interface  showed  an  almost  flat  spectra  while  the  volume  showed  a 
drop-off  with  increaasing  wavenumber. 
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(IPerturbed^/IUnperturbedl2)  Source  depth:  48m,  Receiver  depth:  40m 


Figure  4.35  -  Normalized  Spectrum  of  Volume  Power  Ratio 
Using  Time-Series  Data 
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V.  CONCLUSIONS  &  RECOMMENDATIONS 


A  theoretical  model  and  corresponding  numerical  simulation  for  shallow  water 
reverberation  has  been  developed  for  spatial  and  sound  speed  perturbations  of  the  bottom 
interface  and  volume,  in  the  context  of  the  parabolic  approximation.  The  main  focus  in 
this  thesis  has  been  the  use  of  the  MMPE  model  for  predicting  the  forward  and  back 
propagations,  and  then  developing  post-processing  routines  in  MATLAB  to  compute  the 
reverberation  loss  and  to  perform  signal  analyses  on  the  derived  reverberation. 

The  development  of  statistical  models  to  introduce  interface  and  volume 
perturbations  into  the  MMPE  model  were  also  completed,  noting  that  density  variations 
to  the  volume  have  yet  to  be  included  to  arrive  at  a  more  comprehensive  model.  This  is, 
however,  not  envisaged  to  be  difficult  and  will  be  left  for  follow-on  work.  The  models 
presented  here  were  tested  and  found  to  be  effective  in  modeling  the  perturbations  typical 
of  different  types  of  bottoms  and  sediment  compositions.  An  rms  value  of  1  m  interface 
perturbation  and  a  15  m/s  sound  speed  variation  were  used.  Most  apparent  was  the  effect 
these  perturbations  had  on  a  one-way  propagation  of  the  MMPE  model.  It  showed 
diffused  bottom  penetration  extending  much  farther  in  range  with  less  noticeable  effects 
in  the  water  column. 

Various  analyses  were  conducted  for  both  perturbed  and  unperturbed 
environments  in  order  to  establish  baselines  for  which  to  see  the  effects  of  perturbation. 
This  was  done  in  computing  the  reverberation  losses  across  a  VLA  based  on  a  single 
source  located  at  range,  r  =  0m ,  and  at  48  m  depth.  Arrival  structures  were  clearly 
discernible  due  to  travel  time  differences,  the  shallowest  receivers  showing  later  arrival 
times.  Reverberation  levels  consistently  showed  ‘201og  r’  drop-offs  regardless  of  the 
presence  of  perturbation.  Correlation  across  depth  was  performed  and  showed  that  the 
interface  decorrelated  more  rapidly  than  for  the  volume.  This  was  also  true  of  the 
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scenario  without  perturbations.  Thus,  the  perturbations  were  not  seen  to  affect  correlation 
significantly. 

The  power  spectrum  of  the  reverberation  and  of  the  ratio  of  the  perturbed  to  the 
unperturbed  data  were  computed.  In  the  power  spectral  analysis  of  the  volume 
reverberation,  showed  less  drop-off  for  the  perturbed  data.  However,  both  interface  and 
volume-perturbed  plots  showed  small  and  large-scale  deviations  from  the  unperturbed 
plots,  providing  evidence  of  the  perturbations.  Broadband  analysis  showed  that  repeating 
structures  were  preserved  while  shorter  scale  perturbations  may  have  augmented  the 
spectra.  The  PRSD  showed  a  rather  flat  spectrum  for  the  interface  reverberation  data 
while  volume  reverberation  showed  roll-off  at  higher  wavenumbers.  This  was  consistent 
in  both  the  CW  and  broadband  analyses  but  the  justification  is  not  presently  known. 

In  essence,  we  have  analyzed  the  influence  of  interface  and  volume  perturbations 
on  CW  and  broadband  reverberation  structure  in  shallow  water.  The  impact  of  the 
interface  perturbations  cannot  be  really  concluded  at  this  time  since  the  spectral  analyses 
did  not  provide  any  observable  characteristics.  The  volume  perturbation  spectrum 
appeared  more  amenable  to  extraction  via  signal  processing  and  should  be  explored 
further. 

With  the  conclusion  of  this  thesis,  recommendations  for  future  work  are: 

•  to  generalize  the  Fortran  implementation  of  the  perturbation  models  with  truly 
random  perturbations  by  incorporating  a  ‘seed’; 

•  to  develop  and  implement  density  fluctuations  into  the  sediment  volume  and 
reverberation  theory; 

•  and  to  perform  second  and  third  order  signal  analyses  on  the  reverberation 
data  in  an  attempt  to  extract  information  on  the  perturbation  spectra. 
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